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1.0  INTRODUCTION  AND  SUMMARY 


1,1  INTRODUCTION 

This  report  specifies  the  software  design  of  an  Acquisition  Logic 
Function  (ALF)  for  a monostatic  radar,  terminal  defense,  ballistic 
missile  defense  (BMD)  system.  A separate  report  evaluates  the  algorithms 
specified  in  this  report. This  work  was  performed  by  Systems  Control, 
Inc.  (SCI)  for  the  Massachusetts  Institute  of  Technology  Lincoln  Laboratory 
(MIT-LL)  as  part  of  Purchase  Order  BB-307  under  Prime  Contract 
F19628-76-C-0002  during  the  period  April  1975  through  June  1976,  The 
technical  monitor  for  the  work  was  Dr.  G.  Q.  McDowell. 


1.2  PURPOSE 

MIT  Lincoln  Laboratory,  in  support  of  the  Ballistic  Missile  Defense 
Advanced  Technology  Center  (BMDATC).is  engaged  in  evaluating  bulk  filter 
algorithms  and  acquisition  logic  for  the  terminal  defense  environment. 

The  ultimate  goal  of  the  effort  is  to  design,  Implement,  and  validate  (in 
simulation)  a complete  Acquisition  or  Designation  Function  that  wovild  be 
capable  of  undergoing  field  teats  with  a phased  array  radar.  The  primary 
verification  tool  of  the  Acquisition  Function  will  be  a high-fidelity 
simulation  of  the  radar  environment  as  observed  by  a BMD  system.  This 
simulation  (called  the  D&D  Testbed)  will  exist  on  the  BUD  ARC  computer 
facility. 

The  purpose  of  this  report  is  to  specify  the  software  of  the  ALF 
simulation  to  be  used  on  the  ARC  computer  facility  (CDC  7600)  in  con- 
junction with  the  D4D  Testbed. 


1.3  SPECIFICATION  METHOD 

This  section  describes  the  approach  usad  in  preparing  the  software 
specification.  Although  the  ALF  is  Intended  for  Implementation  on  the 
CDC  7600,  this  report  provides  a computer-independent  specification  of 
Che  ALF  software  design. 
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The  general  form  of  this  report  is  similar  to  the  Milestone  format 
used  by  the  USAF  Space  and  Missile  Systems  Organization.  This  format 
has  proven  itself  to  be  very  useful  for  the  design  specification  of  large 
software  systems,  since  it  clearly  defines  program  structure,  data 
structure,  and  data  flow. 

Flowchart  Conventions 


Figure  1.1  shows  the  flowcharting  conventions  used  throughout  this 
report  for  overview  and  detailed  flowcharts.  Overview  flowcharts  show 
the  logic  flow  and  data  flow  associated  with  the  ALF  and  its  subsystems. 
The  detailed  flowcharts  show  the  logic  flow  for  each  of  the  subroutines 
defined  in  this  report. 

D&D  Testbed/ALF  Integration 

The  ALF  software  is  designed  to  receive  radar  returns  from  the  D&D 
Testbed.  The  origin  of  these  returns  is  denoted  throughout  this  report  as 
the  Matched  Filter  Processing  Function.  The  Testbed  and  the  ALF  will  be 
executed  alternately,  simulating  the  generation  and  processing  of  radar 
returns.  Figure  1.2  shows  an  overview  flowchart  of  this  process.  The 
ALF  software  design  permits  either  open  loop  or  closed  loop  operation  with 
the  Testbed,  as  desired.  In  the  closed  loop  mode  additional  bursts  are 
requested  via  the  Waveform  Request  file  as  shown  in  Figure  1.2.  Through- 
out the  report  this  data  is  described  as  being  sent  to  the  Radar  Scheduler 
Function,  as  would  be  the  case  in  field  implementation.  In  the  open 
loop  mode,  returns  are  generated  in  a prearranged  sequence,  or  recorded 
real  returns  are  used. 

The  Testbed  and  the  ALF  programs  and  associated  data  require  more 
core  than  is  available  in  the  small  core  memory  of  the  CDC  7600.  There- 
fore, some  method  must  be  provided  to  allow  alternate  execution  and  commu- 
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Standard  overlay  processing  (using  LCM  or  disc) 


• Control-point  to  control-point  processing  (each 

program  appearing  as  a separate  job  to  the  operating 
system) 

Development  of  an  integration  approach  requires  a tradeoff  study  involv- 
ing run  time  versus  data  storage  requirements.  The  computer-independent 
software  specification  presented  in  this  report  does  not  assume  a 
specific  integration  approach. 

An  additional  program  called  a Bulk  Filter  Exerciser^ will  be  required 
for  integration.  This  program  will  reformat  data  exchanged  between  the 
Testbed  and  ALF.  If  the  overlay  approach  is  used,  an  overlay  executive 
program  will  also  be  required.  Neither  the  Bulk  Filter  Exerciser  nor 
the  overlay  executive  are  specified  in  this  report.  Their  design  is  con- 
sidered to  be  part  of  the  development  and  integration  phase. 

Data  Structure 


Sections  2 through  9 of  this  report  describe  the  ALF  executive  program 
and  the  subsystems  which  implement  the  component  algorithms.  Each  section 
contains  a subsection  titled  Data  Specification  which  describes  the  arrays 
and  variables  that  comprise  the  input  and  output  data  of  the  program 
or  subsystem.  The  input/output  data  of  each  subsystem  are  listed  in 
groups  with  a common  source  or  destination  to  facilitate  defining  COMMON 
blocks  for  FORTRAN  implementation.  Data  names  use  FORTRAN  Implicit  typing: 
Integer  variables  begin  with  letters  I through  N;  other  variables  are 
real.  For  each  item  or  array,  the  Data  Specification  subsections  provide 
an  engineering  symbol,  which  is  used  in  the  text  and  equations,  and 
a program  symbol,  which  is  used  when  discussing  data  flow,  array  sizes, 
etc. 


A "Parameters"  data  file  is  associated  with  each  algorithm,  e.g., 
the  MARKING  ALGORITHM  PARAMETERS  file.  In  each  case,  such  a file  contains 
parameters  used  only  in  that  algorithm. 


1 

Memory  Management  || 

The  ALF  data  base  (inputs,  outputs,  and  working  files)  is  described  ■ 

in  this  report  as  a set  of  arrays  and  variables.  Because  of  the  • 

large  size  of  the  ALF  data  base,  the  use  of  mass  memory  will  doubtless  . 

be  required  for  some  of  the  data.  When  implementing  the  ALF  on  the  I 

CDC  7600,  data  can  reside  in  small  core  memory  (SCM) , large  core  memory 
(LCM) , or  on  disc.  When  operating  with  the  D&D  Testbed,  data  can  be 
saved  from  one  AIS  execution  to  the  next  by  storing  it  in  LCM  or  on  disc. 

The  location  of  each  array  in  this  memory  heirarchy  is  not  specified  in 
this  report  and  is  considered  to  be  one  of  the  tasks  in  the  development  phase. 
Initial  estimates  for  large  array  sizes  are  provided  in  Section  10.0, 

Memory  Requirements. 

Data  Flow 


Data  flow  is  described  at  three  levels: 

e Between  ALF  and  other  functions 
e Between  ALF  subsystems 
e Between  subroutines  of  a subsystem 

The  description  of  data  flow  is  found  in  the  Data  Specification 
section  for  each  subsystem  and  in  the  overview  flowcharts  and  associated 
text. 

Structured  Programs 

The  ALF  consists  of  a set  of  seven  algorithms  which  are  controlled 
by  the  executive  program.  Each  algorithm  is  implemented  by  a set  of 
subroutines  with  each  subroutine  performing  a single  functional  task. 

The  Individual  subroutines  are  described  in  sufficient  detail  to  identify  j 

their  basic  structure  and  purpose.  Top-down  programming  techniques  are 

used  so  that  logic  flow  proceeds  in  a unidirectional  manner  through  each  | 

subroutine. 

I 
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Coding  Details 


No  attempt  is  made  to  define  each  line  of  code  required  by  each 
subroutine.  Flowcharts  and  text  describe  the  basic  process  performed 
in  each  subroutine.  When  computations  are  to  be  made,  the  flowcharts 
refer  to  equation  numbers  in  the  adjoining  text. 

1.4  EXTERNAL  FUNCTIONS 

This  section  describes  the  other  software  packages  (called  external 
functions)  with  which  the  ALF  is  designed  to  operate. 

The  ALF  simulation  is  designed  to  interact  with  four  external 
functions. 

The  Matched  Filter  Processing  Function  provides  radar  video  and 
monopulse  data  in  a given  beam  for  specific  burst  transmission  and 
threat  conditions.  In  the  simulation,  this  function  would  be  implemented 
by  the  D&D  Testbed.  The  ALF  and  the  D&D  Testbed  would  interact  in  either 
an  open  loop  mode  (predetermined  transmission  sequence)  or  a closed  loop 
mode  (ALF  controls  the  transmission  sequence) . 

The  Radar  Scheduler  Function  controls  the  time  scheduling  for  the 
bursts  in  each  beam.  This  function  contains  an  a priori  search  pattern 
for  a set  of  beams.  In  addition  to  scheduling  the  search  transmissions, 
the  Radar  Scheduler  Function  controls  the  transmission  of  bursts  requested 
by  ALF.  In  the  simulation,  this  function  would  presumably  be  contained 
in  the  D&O  Testbed. 

The  Resource  Allocation  Function  allocates  BMD  system  resources  (e.g., 
burst  transmissions)  among  the  BMD  subsystems  (e.g.,  ALF  and  Track  Func- 
tion) . Future  versions  of  ALF  will  include  a Parameter  Selection  Algori- 
thm for  selecting  among  a set  of  designation  sequences  (see  Section  4.0). 
This  algorithm  will  Interface  with  the  Resource  Allocation  Function. 


The  Track  Function  controls  radar  tracking  o£  the  designated  oOjects 
and  periodically  provides  updated  state  estimates  for  chose  objects.  This 
data  is  used  by  the  Known  Object  Recognition  Algorithm  for  removing  re- 
dundant designations.  If  this  function  is  not  implemented  in  the  simula- 
tion, the  Known  Object  Recognition  Algorithm  will  still  remove  redundant 
objects  by  comparing  designation  files  for  each  beam. 

1.5  OVERVIEW  DISCUSSION 

This  section  provides  a summary  description  of  the  ALF  component 
algorithms. 

The  framework  used  in  designing  the  ALF  is  shown  in  Figure  1.3. 

A 

The  ALF  is  divided  into  seven  component  algorithms  (shown  in  heavy 
outlines).  This  division  is  somewhat  arbitrary;  however,  it  aids  in 
defining  the  requirements  of  the  ALF  and  Che  interfaces  with  the  external 
functions;  radar  scheduler,  resource  allocator,  and  track  function.  The 
arrows  in  Figure  1.3  represent  the  flow  of  information  between  algorithms. 

The  ALF  executive  program  determines  Che  sequence  of  processing 
steps  in  the  ALF  and  passes  data,  as  required,  to  the  various  algorithms. 

The  Waveform  Request  Algorithm  provides  (1)  the  waveform  parameters 
for  each  burst  transmission  and  (2)  the  range*  range-race  acceptance  window 
for  processing  each  beam.  For  the  ALF  design  specified  in  this  report, 
all  of  these  parameters  are  calculated  a priori  and,  hence,  the  Waveform 
Request  Algorithm  simply  accesses  a file  containing  the  information. 


The  ALF  Marking  Algorithm  processes  radar  returns  supplied  by  the 
Matched  Filter  Processing  Function.  This  algorithm  locates  peaks  in 
range  and  range-race  chat  exceed  the  detection  threshold. 

*The  ALF  specified  in  this  document  includes  the  bulk  filter-track 
initiation  algorithm  (l.e.,  the  marking  algorithm,  the  coincidence 
detection  algorithm,  and  the  designation  sequence  combining  algorithm) . 
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The  Coincidence  Detection  Algorithm  processes  the  list  of  marks  from 
two  burst  waveforms  to  remove  range  ambiguities  introduced  by  the  waveform. 
The  SCI  coincidence  detection  algorithm  is  specified  in  this  report. 

The  Coherent  Double  Gating  Algorithm  performs  additional  awbiguity 
removal  and  clutter  rejection;  it  also  provides  initial  state  estimates 
in  range  and  range-rate  for  the  Track  Function.  The  CDG  Algorithm  does  not 
process  angle  measurements. 

The  Parameter  Selection  Logic  executes  the  real-time  selection  of 
ALF  parameters  based  on  the  measured  environment.  The  current  algorithm 
specifies  a predetermined  N-M:N  designation  sequence. 

The  Object  Beam  Position  Logic  computes  the  angular  position  of 
objects  within  the  beam  during  the  designation  sequence. 

The  Known  Object  Recognition  (KOR)  Algorithm  tests  new  designations 
against  targets  in  track  to  prevent  the  initiation  of  redundant  tracks  on 
the  same  target.  The  KOR  algorithm  accomplishes  both  the  "adjacent  beam 
correlation"  and  the  "designation-to-track  correlation"  operations. 


J 


RESOURCE 

ALLOCATION 

FUNCTION 


1.6  MULTI-BEAM  CAPABILITY 


This  section  describes  the  ability  of  the  ALF  to  process  data  from 
a number  of  beams. 

The  ALF  is  designed  to  control  the  designation  process  concurrently 
in  a nvanber  of  beams.  The  transmissions  in  the  beams  are  interleaved  in 
time  in  a manner  determined  by  the  scan  pattern  of  the  Radar  Scheduler 
Function. 

The  multi-beam  capability  is  provided  by  retaining  information  on 
the  objects  in  each  active  beam  from  transmission  to  transmission.  An 
active  beam  is  defined  as  a beam  which  has  entered  the  designation 
sequence. 

The  Known  Object  Recognition  Algorithm  is  a significant  feature  of  the 
multi-beam  capability.  This  algorithm  removes  redundant  designations  in 
adjacent  beams,  i.e.,  designations  on  the  same  object  in  two  or  more 
beams . 
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2.0  ALF  EXECUTIVE  PROGRAM 


2 . 1 PURPOSE 

The  ALF  executive  program  controls  the  execution  of  the  ALF  soft- 
ware. The  executive  calls  each  of  the  subroutines  which  implement  ALF 
subsystems. 

2.2  GENERAL  DESCRIPTION 

This  section  describes  the  sequence  in  which  the  ALF  algorithms 
are  executed  and  the  data  flow  between  algorithms. 

The  structure  of  the  ALF  and  the  basic  purpose  of  each  of  its  sub- 
systems is  described  in  Section  1. 

In  each  pass  through  the  ALF  software,  radar  data  is  processed  for 
one  burst  from  one  beam.  Various  data  files  are  retained  concerning  the 
status  of  that  beam;  this  data  is  used  and  updated  when  Che  next  burst 
from  chat  beam  is  processed. 

The  ALF  is  called  by  an  executive  program  whenever  radar  data  is 
available  for  a burst  from  a beam.  The  Matched  Filter  Processing  Function 
provides  Che  radar  video  and  monopulse  data  for  the  beam.  The  beam  number 
and  ocher  data  which  identify  the  source  and  nature  of  Che  radar  returns 
are  also  provided.  If  this  data  exists  on  disc,  the  ALF  executive  reads 
this  data  and  scores  it  in  core  memory  (probably  LCM) . 

The  ALF  design  assumes  that  Che  Radar  Scheduler  Function  will  initially 
cause  Che  Matched  Filter  Processing  Function  to  provide  radar  and  monopulse 
data  on  a predetermined  schedule  for  each  beam.  This  process  simulates  Che 
transmission  of  search  bursts.  The  purpose  of  a search  burst  is  to  determine 
whether  there  arc  any  objects  in  a beam  and,  if  so.  to  provide  an  initial 
sec  of  marks  for  that  beam. 

Figure  2.1  is  an  overview  flowchart  for  the  ALF  software. 
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If  Che  first  burst  of  a pair  is  being  processed,  Figure  2.1  shows  Che 
Marking  Algorithm  (MRKALG)  is  called  to  process  the  BURST  1 VIDEO  file  and 
to  pick  peaks  in  range  and  range-rate.  The  location  of  these  marks  is  stored 
in  Che  BURST  1 MARKS  file.  If  there  are  any  marks  and  if  Che  burst  was  a 
search  burst,  Chen  a second  burst  is  requested  by  the  Waveform  Request 
Algorithm  2-3  msec  after  the  search  burst.  If  there  are  no  marks  on  a search 
burst,  no  additional  burst  is  requested,  and  the  next  transmission  in  that 
beam  will  be  determined  by  Che  a priori  schedule  of  the  Radar  Scheduler  Function. 


I 

I 
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If  the  ALF  is  processing  the  second  burst  of  a pair,  processing 
follows  Che  logic  path  shown  on  Che  right  hand  side  of  the  first  page  of 
Figure  2.1.  Marking  is  performed  on  the  BURST  2 VIDEO  file  and  the 
marks  are  scored  in  the  BURST  2 MARKS  file.  The  Coincidence  Detection 
Algorithm  then  examines  marks  and  video  from  both  bursts  to  produce  a list 
of  COINCIDENCE  DETECTION  MARKS. 

If  the  burst  pair  number  equals  or  exceeds  a parameter  M (defined 
below)  the  Object  Beam  Position  Algorithm  (ANGALG)  is  called  to  compute 
angular  positions  of  all  objects  whose  marks  passed  coincidence  detec- 
tion. This  data  is  stored  with  each  mark  in  the  COINCIDENCE  DETECTION 
MARKS  file. 

If  the  burst  pair  being  processed  is  the  first  pair  of  a designation 
sequence,  a test  is  made  to  determine  if  any  marks^  passed  coincidence 
detection.  If  no  marks  passed,  processing  is  terminated  and  the  beam 
reverts  to  the  search  mode.  If  there  are  coincidence  detection  marks, 
the  Parameter  Selection  Algorithm  is  called  to  determine  an  N-M:N  desig- 

it 

nation  sequence  for  the  beam.  The  Coherent  Double  Gating  Algorithm 
(CDGALG)  is  called  and  initialized  (for  the  set  of  coincidence  detection 
marks)  to  form  the  CDG  TRACK  FILE.  The  Waveform  Request  Algorithm  is  then 
called  to  request  the  next  burst  pair  50  msec  later.  Once  the  designa- 
tion sequence  of  N burst  pairs  has  begun,  burst  pairs  are  requested  rather 
than  single  bursts. 


In  the  CDG  algorithm  a strategy  which  requires  n-M  detections  in  a total 
of  N burst  pair  transmissions  is  denoted  by  N-M:N. 


I 
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video  data  from  successive  cransmlssloas  are  processed  as  described 
above  up  to  the  point  marked  2A  on  Figure  2.1.  At  this  point,  the 
Coherent  Double  Gating  Algorithm  Is  called  to  update  the  object  tracks 
from  the  previous  transmission  (Page  2 of  Fig.  2.1).  Those  objects 
whose  tracks  pass  the  gating  test  are  retained  (In  updated  form)  In  the 
CDG  TRACK  FILE.  Then  If  there  are  tracks  remaining  in  the  beam  and  the 
designation  sequence  Is  not  completed,  another  burst  pair  Is  requested 
by  the  Waveform  Request  Algorithm.  If  there  are  no  tracks  remaining, 
the  beam  reverts  to  the  search  mode.  If  there  are  tracks  remaining  at 
the  end  of  the  designation  sequence  (after  burst  N) , the  Known  Object 
Recognition  Algorithm  (KORALG)  Is  called  to  remove  redundant  designations 
within  the  beam.  The  designations  corresponding  to  the  adjacent  beams  in 
DESIGNATIONS  file  are  examined  to  remove  any  designations  from  the  current  beam 
which  are  redundant  between  beams.  The  KOR  Algorithm  also  removes  designations 
from  the  current  beam  which  are  Indistinguishable  from  objects  already  In 
track.  Those  objects  within  the  beam  which  pass  the  KOR  test  are  passed 
to  the  Track  Function  in  the  DESIGNATIONS  File. 

Figure  2.2  shows  the  sequence  of  events  for  processing  a single  beam. 

Each  horizontal  sequence  describes  the  generation  and  processing  of  returns 
for  one  burst  in  that  beam.  The  sequences  in  Figure  2.2  are  shown  in  the 
time  order  of  occurrence  (from  the  top  down).  These  sequences  depict  each 
major  logic  path  through  the  overview  flowchart  in  Figure  2.1. 

In  Figure  2.2  a reference  to  burst  (k,  1)  means  the  burst  (first 
or  second)  of  burst  pair  k,  where 

1 < k <.  N. 

Burst  (1,1)  is  the  search  burst  and  la  requested  by  the  Radar  Scheduler 
Function.  The  presence  of  marks  on  a search  burst  initiates  the  acquisi- 
tion sequence,  i.e.,  burst  (1,2)  is  requested,  as  shown  in  the  second 
sequence  in  Figure  2.2. 
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Figure  2.1  ALF  Software  Overview  Flowchart 
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Figure  2.1  Software  Overview  Flowchart  (Continued) 
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Burst  transmissions  for  various  beams  are  interleaved  in  time. 
Therefore  the  single  beam  processing  shown  in  Figure  2.2  is  interleaved 
with  the  processing  of  other  beams,  since  bursts  are  processed  in  the 
order  in  which  they  are  received. 

The  BEAM  STATUS  file  contains  information  which  describes  the 
current  status  of  each  beam  which  is  either  undergoing  or  has  completed 
designation  processing.  The  file  entry  for  a beam  is  updated  by  the 
ALF  executve  program  each  time  a burst  is  processed  for  that  beam. 
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2.3  DATA  SPECIFICATION 

I 

This  section  discusses  the  data  interchanged  between  the  ALF 
executive  and  the  external  functions.  A bookkeeping  File  (BEAM  STATUS)  | 

is  also  discussed. 

I 

2.3.1  Input  Data 

The  following  data  is  input  from  the  Matched  Filter  Processing 
Function: 


Program 

Symbol 

Engineering 

Symbol 

Description 

lU 

- 

Index  for  beam  U coordinate 

IV 

- 

Index  for  beam  V coordinate 

NBPAIR 

- 

Burst  pair  number  for  current  burst 

NBURST 

- 

Location  of  the  current  burst  within 
a burst  pair  (NBURST  = 1,2) 

RS 

^S 

Range  corresponding  to  start  of 
video  data  processing  (m) 

RE 

Range  corresponding  to  end  of  video 
data  processing  (m) 

DELR 

A 

Range  sample  spacing  (m) 

DS 

°S 

Range-rate  corresponding  to  lowest 
Doppler  channel  for  which  video  data 
is  processed  (m/s) 

DE 

°E 

Range-rate  corresponding  to  highest 
Doppler  channel  for  which  video  data 
is  processed  (m/s) 

DEED 

Range-rate  sample  spacing  (Doppler 
channel  spacing)  (m/s) 

VIDEO(K,I) 

r(k;l) 

Amplitude  of  video  at  range  Rg  + (k-l)A 
and  range-rate  Dg  + (l-l)c.  Note  that 
k“l  at  Rg  and  i*l  at  D^. 

UBSTl(K.I) 

Au(k, i) 

Monopulse  measurement  of  the  U coordinate 
indexed  like  r(k,i) 

VBSTKK.I) 

Av(k,i) 

Monopulse  measurement  of  the  V coordinate 
Indexed  like  r(k,i) 
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The  following  data  is  input  from  the  TRACK  FUNCTION  via  the 


OBJECT  TRACK  file: 


Program 

Symbol 


Engineering 
Symbol 


Description 


NTRK 


Number  of  objects  in  the  array  TRACKS 


TRACKS (NTRK, I)  - List  of  the  most  recent  data  on  the 

objects  in  track 

i * 1 Range  estimate  (m) 

i • 2 Range-rate  estimate  (m/sec) 

2 

1 ■ 3 Variance  of  range  estimate  (m  ) 

i * 4 Covariance  of  range  and2 

range-rate  estimates  (m  /sec) 

i ■ 5 Variance  of  the  range-rate 
estimate  (m2/sec2) 

i • 6 Angular  coordinate  U 

i • 7 Angular  coordinate  V 

i ■ 8 Time  of  estimation  (msec) 


The  following  data  is  input  via  the  BEAM  STATUS  File: 


ICUR 

- 

Index  to  array  entries  for  current  beam 

IDES(ICUR) 

- 

Flag  indicating  designation  status 
of  beam 

IDES (ICUR)  ■ 1 Designation  complete 
IDES(ICUR)  » 0 Designation  incomnlete 

TBURST(ICUR) 

- 

Time  of  last  burst  transmission  in 
the  beam  (msec) 

lU(ICUR) 

- 

Index  for  beam  U coordinate 

TV  (ICUR) 

- 

Index  for  beam  V coordinate 

NBPAIR(ICUR) 

- 

Number  of  burst  pair  being  processed 
in  the  beam 

NBURST(IOTR) 

- 

Number  of  burst  being  processed 
in  the  beam 

NBEAM 

- 

Number  of  beams  currently  listed 
in  the  BEAM  STATUS  file 

2.3.2  Output  Data 


The  following  data  is  output  to  the  Radar  Scheduler  Function: 


Program 

Symbol 

NBREQ 

TNEXT 


TSEP 


TIPSl 


Engineering 
Symbol 


Description 

Number  of  bursts  requested  (1  or  2) 

Transmission  time  for  next  requested 
burst  in  the  beam  (typically  50  msec 
after  current  burst  pair) 

Time  separation  between  bursts  in 
the  requested  burst  pair 

Interpulse  spacing  in  the  next 
requested  burst  (psec) 


I 

I 


TIPS2 

lU 

IV 


Interpulse  spacing  in  the  second  burst 
of  requested  burst  pair  (usee) 

Index  for  beam  U coordinate 

Index  for  beam  V coordinate 


The  following  data  is  output  to  the  Track  Function  via  the 
DESIGNATIONS  file: 


NDES  - Number  of  objects  in  the  array  DESIG 

DESIG(NDES , I)  - List  of  objects  designated  in  Che 

beams  already  processed  in  the 
current  scan  of  the  group  of  beams. 


The  updated  BEAM  STATUS  file 


1 ■ 1 Range  estimate  (m) 

2 

i ■ 2 Range-rate  estimate  (m  ) 

2 

i ■ 3 Variance  of  range  estimate  (m  ) 

i ■ 4 Covariance  of  range  and, 

range-rate  estimates  (m*'/sec) 

i - 5 Variance  of  range-rate 
estimate  (m2/sec2) 

1 ■ 6 Angular  coordinate  U 

i ••  7 Angular  coordinate  V 

1 ■ 8 U index  for  beam  of  designation 

1 ■ 9 V index  for  beam  of  designation 

is  output  by  the  ALF  executive. 
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2.4  MATHEMATICAL  RELATIONSHIPS 


The  mathematical  relationships  embodied  In  the  ALF  algorithms  are 
described  In  the  sections  which  specify  the  Individual  algorithms.  The 
executive  program  provides  the  logic  for  sequential  execution  of  the 
algorithms. 
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3.0  WAVEFORM  REQUEST  ALGORITHM 


3.1  PURPOSE 

The  Waveform  Request  Algorithm  requests  the  Radar  Scheduler 
Function  to  transmit  specific  waveforms  required  for  bulk  filtering. 

It  also  provides  parameters  required  by  other  algorithms  for  process- 
ing the  radar  return  signal. 

The  WAVEFORM  REQUEST  file  is  generated  only  when  the  D&D  Testbed/ 

ALF  is  operating  in  the  closed  loop  mode,  i.e.,  only  when  the  ALF  determines 
what  burst  will  be  transmitted  next. 

3.2  GENERAL  DESCRIPTION 

Figure  3.1  is  the  overview  flowchart  for  the  Waveform  Request  Algorithm 
(WAVALG) . The  BEAM  STATUS  file  Indicates  which  beam  has  just  been  processed, 
the  transmission  time  of  the  received  burst,  and  the  burst  pair  and  burst 
number.  From  this  data,  the  Waveform  Request  Algorithm  computes  the  trans- 
mission time  and  interpulse  spacing  for  the  next  burst  or  burst  pair  in  the 
beam.  These  computed  burst  parameters  are  stored  in  the  WAVEFORM  REQUEST  file 
(along  with  the  beam  indices,  lUREQ  and  IVREQ)  for  use  by  the  Radar  Scheduler 
Function. 

The  Waveform  Request  Algorithm  chooses  a different  Interpulse  spacing 
from  a prestored  table  for  each  successive  requested  burst.  The  pulse 
spacing  is  varied  to  provide  additional  range  ambiguity  removal  in  the 
Coherent  Double  Gating  Algorithm. 
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3.3  DATA  SPECIFICATION 


3.3.1  Input  Data 


The  following  data  is  input  from  the  ALF  Executive  via  the  BEAM  STATUS 
File; 


Program 

Symbol 

Engineering 

Symbol 

Description 

ICUR 

- 

Index  to  array  entries  for  current  beam 

IDES(ICUR) 

Flag  indicating  designation  status  of  beam 
IDES(ICUR)  ■ 1 Designation  complete 
IDES(ICUR)  ■ 0 Designation  incomplete 

TBURSTdCUR) 

- 

Time  of  last  burst  transmission  in  the 
beam  (msec) 

lu(iom) 

- 

Index  for  beam  U coordinate 

IV(ICUR) 

- 

Index  for  beam  V coordinte 

NBPAIR(iaTl) 

- 

Number  of  burst  pair  being  processed  in 
the  beam 

NBURST(ICUR) 

- 

Number  of  burst  being  processed  in  the  beam 

NBEAM 

- 

Number  of  beams  currently  listed  in  the 
BEAM  STATUS  file. 

The  following  data  is  input 

via  the  WAVEFORM  PARAMETERS  file: 

VALIPS(I.J) 

- 

Prestored  values  of  interpulse  spacing 
(usee) 

I Burst  pair  number 

J Burst  number 

TB 

^B 

Time  between  bursts  (msec) 

TP 

^P 

Time  between  burst  pairs  (msec) 

BEAMS  (lU,  IV. L) 

- 

List  of  range,  range-rate  processing  limits 
for  each  beam  in  the  search  sector 

lU  Beam  U index 
IV  Beam  V index 
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R 


SW 


L ■ 1 Range  corresponding  to  closest 
edge  of  acquisition  band  in  the 
beam  (m) 

L ” 2 Range  corresponding  to  the 

furthest  edge  of  acquisition 
band  in  the  beam  (m) 

L “ 3 Range-rate  corresponding  to  the 
lowest  Doppler  channel  to  be 
processed  for  the  beam  (m/sec) 

L - 4 Range-rate  corresponding  to  the 
highest  Doppler  channel  to  be 
processed  for  the  beam  (m/sec) 


3.3.2  Output  Data 


The  following  data  is  output  to  the  Radar  Scheduler  Function  via 
the  WAVEFORM  REQUEST  File: 


Program 

Symbol 

NBREQ 

TNEXT 

TSEP 

flPSl 

TIPS  2 

lUREQ 

IVREQ 


Engineering 

Symbol Description 

Number  of  bursts  requested  (1  or  2) 

Transmission  time  for  next  requested  burst 
in  the  beam  (msec) 

Tg  Time  separation  between  bursts  in  the 

requested  burst  pair  (msec) 

- Interpulse  spacing  in  the  next  requested 
burst  (usee) 

- Interpulse  spacing  in  the  second  burst  of 
requested  burst  pair  (ysec) 

Index  for  beam  U coordinate 

- Index  for  beam  V coordinate 


w 
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3.3.3  Parameters 


Table  3.1  shows  recommended  values  of  interpulse  spacing  for 
three  successive  burst  pairs  in  a beam,  where  the  nominal  pulse  spacing 
is  5.0  usee.  Data  is  stored  in  array  VALIPS(I.J).  By  varying  the  pulse 


Table  3.1 

RECOMMENDED  INTERPULSE  SPACING  VALUES  (nominal-5 . 0 usee) 


BURST  PAIR 
SEQUENCE  NUMBER 

BURST  NUMBER  WITHIN 
A BURST  PAIR 

1 2 

1 

5.0  sec 

4.7  sec 

2 

5.2 

4.9 

3 

5.1 



4.8 

spacing  on  successive  burst  pairs,  the  Coherent  Double  Gating  Algorithm 
can  reject  ambiguities  which  have  passed  through  the  Coincidence  Detec- 
tion Algorithm.  For  designation  sequences  longer  than  3 burst  pairs, 
the  sequence  of  interpulse  spacings  in  Table  3.1  is  repeated  as  many  times 
as  necessary. 

Other  parameters  are  as  follows: 

T Time  between  bursts  (nominally  2 to  3 msec) 

D 

Tp  Time  between  burst  pairs  (nominally  50  msec) 

The  two  parameters  Tg  and  Tp  are  defined  as  shown  in 
Figure  3.2. 
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Figure  3.2  Burst  Transmission  Timing  Definitions 
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Figure  3-3  Waveform  Request  Algorithm  Detailed  Flowchart 
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4.0  PARAMETER  SELECTION  ALGORITHM 


I 

^ 4.1  PURPOSE 

The  Parameter  Selection  Algorithm  executes  the  real-time  selection  of 
ALF  parameter  values  based  on  the  measured  environment.  Since  the  develop- 
ment of  an  algorithm  for  choosing  the  N-M:N  designation  sequence  has  not 
I been  completed,  the  current  Parameter  Selection  Algorithm  specifies  a pre- 

determined N-M;N  sequence. 

Future  versions  will  include  an  algorithm  for  selecting  among  a set 
I of  designation  sequences. 

I 4.2  GENERAL  DESCRIPTION 

. Figure  4.1  is  the  overview  flowchart  for  the  Parameter  Selection 

I Algorithm  (PARALG) . The  Parameter  Selection  Algorithm  contains 

the  predetermined  values  of  M and  N for  the  designation  sequence.  These 
I values  are  passed  to  the  CDG  Algorithm  for  each  beam  via  the  CDG  PARAMETERS  file. 

I 4.3  DATA  SPECIFICATION 

I 4.3.1  Input  Data 

I The  following  data  is  input  to  the  Parameter  Selection  Algorithm: 


Program 

Engineering 

Symbol 

Description 

lU 

- 

Index 

for  beam  U coordinate 

IV 

Index 

for  beam  V coordinate 
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4.3.2  Output  Daca 


The  following  data  is  output  to  the  CDG  Algorithm  via  the  CDG  PARAMETERS 

file: 


Engineering 

Symbol  Description 

N Number  of  transmissions  in  designation 

sequence . 

M Number  of  misses  allowed  in  designation 

sequence. 

4.3.3  Parameter  Settings 

Designation  sequences  with  M = 1,  N =■  6;  M ” 2,  N * 6;  or  M = 1, 

N ■ 5 are  recommended. 


Program 

Symbol 

NTRAN 

NMISS 


I 
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5.0  MARKING  ALGORITHM 


5.1  PURPOSE 

The  Marking  Algorithm  operates  on  the  sampled  radar  video  from  several 
Doppler  filters  to  generate  marks,  i.e.,  estimates  of  the  two~element  state 
of  range  and  range-rate  for  potential  objects.  Basically,  the  Marking 
Algorithm  selects  peaks  in  the  radar  video  response  (over  a portion  of  the 
range,  range-rate  plane)  which  have  amplitudes  characteristic  of  potential 
RVs.  For  each  of  these  peaks  the  algorithm  estimates  the  location  in  the 
range,  range-rate  plane  and  saves  it  for  further  processing. 


5.2  GENERAL  DESCRIPTION 

The  range,  range-rate  Marking  Algorithm  described  in  this  report  is 
designed  to  operate  on  sampled  video  (sampling  rate  in  range  equal  to 
twice  the  nominal  bandwidth)  from  Doppler  filters  spaced  at  305  m/s  in 
range  rate.  The  algorithm  should  operate  satisfactorily  with  Doppler 
filter  spacing  as  large  as  one-half  of  the  Doppler  resolution.  These 
restrictions  on  sampling  races  in  range  and  Doppler  were  imposed  to 
eliminate  the  requirement  for  sophisticated  interpolation  techniques. 

Since  the  input  data  are  assumed  to  be  calibrated  in  range  and  range- 
race,  all  opera.:ions  are  described  in  these  units  (rather  than  time  and 
Doppler  frequency).  Video  samples  are  provided  for  all  Doppler  filters 
within  the  RV  range-rate  acceptance  region  plus  one  filter  on  each  side 
of  the  region.  This  allows  marks  to  be  generated  in  one  filter  beyond 
each  edge  of  the  RV  range-rate  acceptance  region.  The  acceptance  of  marks 
in  one  filter  beyond  the  RV  range-rate  acceptance  region  reduces  the 
probability  of  missing  an  RV  due  to  the  Doppler  resolution  of  the  waveform. 
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A range,  range-rate  mark  is  generated  by  either  of  the  following 
two  criteria: 

1.  A peak  in  range-rate  occurs  at  the  same  sample  location 

as  a peak  in  range,  and  the  video  amplitude  at  that  sample 
exceeds  the  detection  threshold  (Type  1 mark) ; 
or 

2.  A peak  in  range-rate  occurs  one  range- rate  sample  away 
from  a peak  in  range,  the  video  amplitude  at  both  samples 
exceeds  the  detection  threshold,  and  a mark  generated  by 
the  first  criterion  does  not  lie  within  plus  or  minus  one 
sample  spacing  in  range  and  plus  or  minus  one  sample  spacing 
in  range-rate  (Type  2 mark) . 


The  second  criterion  for  generating  a mark  reduces  leakage  caused  by 
interference  from  fragments  or  fragment  ambiguities  at  lower  range-rates. 
The  final  estimates  of  range  and  range-rate  for  each  mark  are  computed  by 
linear  interpolation  of  the  differences  between  sample  amplitudes. 

Figure  5.1  is  an  overview  flow  chart  showing  sequence  of  events  and 
data  flow  for  the  marking  algorithm.  The  range,  range-rate  marking 
algorithm  first  examines  video  data  to  locate  peaks  in  range  which  exceed 
the  detection  threshold,  and  then  proceeds  to  locate  peaks  in  range-rate 
at  each  of  these  range  locations.  The  resulting  range,  range-rate  marks 
are  tested  for  satisfaction  of  the  marking  criteria. 
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I 


MRKALG 


i 


I 


I 


Figure  5.1  Marking  Subsystem  Overview  Flowchart 


I 
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5.3  DATA  SPECIFICATION 


5.3,1  Input  Data 


The  following  data  is  input  via  the  MARKING  PARAMETERS  file: 


Program 

Engineering 

Symbol 

Symbol 

Description 

RSW 

^SW 

Range  corresponding  to  closest  edge  of 
acquisition  band  in  the  beam  (m) 

REW 

Range  corresponding  to  furthest  edge  of 
acquisition  band  in  the  beam  (m) 

RREF 

"o 

Reference  range  for  specification  of 
detection  threshold  (ra) 

THRESH 

"o 

Detection  threshold  (specified  as  a radar 
cross  section)  at  the  reference  range 

RQ(m2) . 

The  following  parameters 
BURST  2 VIDEO  file: 

are  input  via  either  the  BURST  1 VIDEO  or 

RS 

"s 

Range  corresponding  to  start  of  video 
data  processing  (m) 

RE 

Range  corresponding  to  end  of  video 
data  processing  (m) 

DELR 

A 

Range  sample  spacing  (m) 

DS 

°s 

1 

Range-rate  corresponding , to  lowest  Dopplei 
channel  for  which  video  data  is 
processed  (m/s) 

DE 

°E 

Range-rate  corresponding  to  highest 
Doppler  channel  for  which  video  data 
is  processed  (m/s) 

DEED 

C 

Range-rate  sample  spacing  (Doppler  channel 
spacing)  (m/s) 

VIDEO(K,I)  r(lt;i)  Amplitude  of  video  at  range  Rg  + (k  -1)A 

and  range-rate  Dg  + (i  - l)c.  Note  that 

k ■ 1 at  Rg  and  i ■ 1 at  Dg  . 
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5.3.2  Output  Data 


The  following  data  is  output  via  either  the  BURST  1 MARKS  or 
the  BURST  2 MARKS  file: 


Program 

Symbol 

NM 

RRDOT(Q.M) 


Engineering 
Symbol 


RR(q,m) 


Description 

Number  of  range,  range-rate 
marks  in  array  RRDOT 

Range-ordered  list,  indexed  by  q, 
of  range,  range-rate  marks  generated 
for  the  current  burst 

m - 1 Range  (m) 

m ” 2 Range-rate  (m/sec) 

m ■ 3 Marking  flag 


5.3.3  Parameter  Settings 

The  only  parameters  which  must  be  specified  in  the  range,  range-rate 

marking  algorithm  are  the  detection  threshold,  T^,  and  the  reference  range, 

Rq.  Normally,  the  middle  of  the  range  corresponding  to  the  acquisition 

band  is  used,  i.e.,  R_  ■ 1/2(R„,  + R_.) . However,  since  the  interference 

U bw  ciw 

often  cannot  be  treated  as  white  gaussian  noise,  standard  operating 

characteristics  cannot  be  used  exactly  to  achieve  a desired  P^.  A 

sensitivity  analysis  via  simulation  is  the  most  reliable  way  to  select 

-4  2 

the  threshold.  Values  of  detection  threshold  between  5*10  m (-33  dBsm) 
-4  2 

and  3.16*10  m (-35  dBsm)  have  been  simulated. 
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5.4  MATHEMATICAL  RELATIONSHIPS 


Range  marking  is  performed  for  one  range  sample  at  a time,  starting 
with  range  R^.  At  each  range  sample,  range  marking  is  performed  for  each 
of  the  Doppler  channels.  This  is  accomplished  by  the  following  sequence 
of  operations.  At  range  sample  k,  the  successive  amplitude  differences 
(for  Doppler  channel  i) 

d(k;  i)  • r(k;  i)  - r(k-l;  i)  , 


+ 2 c k < 


+ 1 


(5.1) 


are  computed.  If  d(k;  i)  > 0 and  d(k  + 1;  i)  £ 0 and  r(k;  i)  ^ T^^  , 
a range  mark  is  generated  at  range 


R(j;  i)  - 


(k  - + 


d(k;i) 


d(k;  i)  - d(k+l;  i) 


(5.2) 


where  j indicates  the  order,  in  range,  of  the  mark  for  Doppler  channel  i. 
As  indicated  in  equation  (5.1),  the  marking  is  performed  only  over  the 
range  interval  [R_„,  R^]  • The  detection  threshold  T , is  varied  with 

i!*W  iC 

range  according  to 


T 


k 


(5.3) 


Next,  for  each  range  mark  R(j ; 1),  a check  is  made  to  determine  if 
a peak  in  range-rate  coincides  with  it.  Each  of  the  range  marks  is  checked 
in  order  of  range,  regardless  of  Doppler  channel.  At  the  range  sample 
nearest  to  R(J ; i) , that  is 


k ■ Integer 


R(j;  1)  - Rg 

A 


I 


(5.4) 
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Che  ampllcude  differences 


e(k;  1)  - r(k;  i)  - r(k;  i - 1)  , 
i - 2,  3,  .... 


(5.5) 


are  computed,  vhere 


If  e(k;  p)  > 0 and  e(k;  p+1)  _<  0 and  r(k;  p)  ^ T^,  a range-rate 
mark  is  generated  at  range-rate 


R(k;  i)  - + 


(P  - f); 


P) 


e(k;  p)  - e(k;  p+1) 


(5.6) 


vhere  I indicates  the  order,  in  range-rate,  of  range-rate  peaks  at 
range  sample  k.  If  no  peaks  in  range-race  are  found,  the  range  mark 
R( j ; i)  is  dropped.  If  one  or  more  range-rate  marks  are  generated  at 
range  sample  k,  each  is  tested  to  determine  if  it  coincides  with  the  range 
mark  R(j ; 1) . If 

|R(k;  Z)  - Dg  - (i  - Del  < e , (5.7) 

then  the  range  estimate  R(j;  i)  and  the  range-rate  estimate  R(k;  1) 
are  scored  in  positions  m ■ 1 and  m ■ 2 respectively  of  the  array 
RR(q,  m) , where  q indicates  the  order  in  which  the  range,  range-rate 

e 

marks  are  generated.  Position  m • 3 in  RR(q,  m)  is  sec  as  follows 

, if  lR(k;  D - Dg  - (1  - Del  < I . 

, otherwise 


RR(q,3) 
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After  the  range,  range-rate  marks  RR(q,  m)  are  generated  as  described 
above,  they  are  reordered  by  increasing  range.  Marks  are  then  retained  in 
RR(q,  m)  if  one  of  the  following  two  conditions  is  satisfied: 


5.5  DETAILED  DESCRIPTIONS 


5.5.1  Marking  Algorithm  Executive  Subroutine  (MRKALG) 

5. 5. 1.1  Purpose 

The  executive  subroutine  MRKALG  controls  the  process  of  marking 
video  data  from  one  burst. 

5. 5. 1.2  Detailed  Description 

The  operation  of  the  executive  is  very  straightforward,  as  shown  In 
the  detailed  flowchart  of  Figure  5.2.  The  sort  operation  simplifies  the 
search  for  redundant  marks . 
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MRKALG 


MARK \ 

PRODUCE  LIST  OF 
RANGE,  RANGE-RATE 
MARKS  BASED  ON  MARKING  RULES, 


Figure  5.2  Detailed  Flowchart  for  Marking  Algorithm 
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5.5.2  MARK 


\ 


I 

I 

I 


5. 5. 2.1  Purpose 


MARK  controls  the  process  of  locating  range  and  range-rate  peaks. 
It  then  tests  for  satisfaction  of  the  marking  criterion  (except  for 
redundant  mark  test) . 

5. 5. 2. 2 Input/Output  Data 


Input  Data 


Program  Engineering 

Symbol  Symbol 


Description 


VIDEO(K,I)  r(k,i) 


Amplitude  of  video  at  range  R^  + (K-l)A 
and  range-rate  D^  + (i-Dc 


Output  Data 

RRDOT(Q,M)  RR(q,m)  Range-ordered  list  of  marks 


5. 5. 2. 3 Detailed  Description 

MARK  produces  range,  range-rate  marks  by  applying  the  marking  rules 
specified  in  Section  5.4.  The  detailed  flowchart  of  MARK  is  shown  in 
Figure  5.4.  The  video  data  is  examined  one  range  bin  at  a time;  subroutine 
PEAK  is  called  to  detect  range  peaks.  Subroutine  VPEAK  is  called  to  locate 
range-rate  peaks.  A search  is  then  made  to  match  each  range-rate  mark 
with  the  closest  range  mark.  If  both  marks  occur  in  the  same  range-rate 
bin  (type  1 mark)  a special  flag  is  set  for  that  range,  range-rate  mark. 

The  only  other  marks  retained  are  those  where  range  and  range-rate  marks 
occur  in  adjacent  range-rate  bins  (type  2 mark). 
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I 

Figure  5.3  shows  an  example  of  the  marking  process  for  the  range  bin 
corresponding  to  range  R^.  Range  marks  and  r-nge-rate  marks  are  shown  | 

for  several  samples  where  there  are  threshold  crossings.  The  range, 
range-rate  mark  at  A will  be  flagged  as  a type  1 mark;  B is  a type  2 mark.  j 

The  detailed  flowchart  of  subroutine  MARK  is  shown  in  Figure  5.4.  I 
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R.  RANGE  BIN 
1 


[]]  « Threshold  Crossing 
X • Range  mark 

0 “ Range-rate  mark 

A ” Range,  range-rate  mark 


Figure  ^.3  Example  of  the  Marking  Process 
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Figure  5.4  Subroutine  MARK  Detailed  Flwochart  (Cont'd) 
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5.5.3  PEAK 


5. 5.3.1  Purpose 

PEAK  processes  video  amplitude  data  to  test  for  a peak  in  range  for 
a specific  range  sample  and  doppler  channel. 

5.5. 3. 2 Essential  Input/Output  Data 
Input  Data 

k - range  sample  number 


r(q,i) 


doppler  channel  number 
- video  data  (k-1  £ q £ k+1)  (array  VIDEO) 
detection  threshold  at  range  sample  k 


Output  Data 


PKFLG 


interpolated  range  value  for  range  peak 
at  range  sample  k 

peak  flag 


5. 5. 3. 3 Detailed  Description 

PEAK  first  tests  whether  the  video  sample  exceeds  the  detection 
threshold  for  this  range.  If  so,  successive  amplitude  differences 

are  computed  and  the  peak  test  is  made.  If  a peak  is  found,  its  range 
is  computed  by  interpolation.  Figure  5.5  shows  an  example  of  the 
detection  of  a peak  between  range  samples  and  iri  doppler 

channel  i . 

The  detailed  flowchart  of  subroutine  PEAK  is  shown  in  Figure  5.6. 


. • • 
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PEAK 


Sutroucine  PEAK  Detailed  Flowchart 


Figure  5.6 
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5.5.4  VPEAK 


5. 5. 4.1  Purpose 

VPEAK  processes  video  amplitude  data  to  test  for  peaks  in  range-rate 
at  a specified  range  sample. 

5. 5. 4. 2 Essential  Input/Output  Data 


• Input  Data 
k 

r(k,i) 


range  sample  number 
- video  (array  VIDEO) 


• Output  Data 
R(i) 

NPKS 


interpolated  range- rates  for  range-rate 
peaks  at  this  range  sample  (array  VELMRK) 

number  of  peaks  found 


5. 5. 4. 3  Detailed  Description 


Subroutine  VPEAK  performs  essentially  the  same  peak  detecting 
functions  in  range-rate  as  are  performed  by  PEAK  for  range.  In  addition, 
the  interpolated  range-rate  is  tested  to  determine  whether  it  lies  within 
the  region  of  interest;  if  it  lies  outside  the  region,  the  mark  is  deleted. 

The  detailed  flowchart  of  VPEAK  is  shown  in  Figure  5.7. 
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Figure  5.7  Subroutine  VPEAK  Detailed  Flowchart 
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5.5.5  SORT 


5. 5. 5.1  Purpose 


The  subroutine  SORT  sorts  the  RR  marks  on  range  (ascending  order). 
The  operation  of  this  subroutine  is  self-explanatory. 
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5.5.6  REMOVE 


5. 5. 6.1  Purpose 

REMOVE  searches  the  range-ordered  list  of  marks  and  removes  type  2 
marks  which  are  adjacent  to  type  1 marks,  as  described  in  the  marking  rules. 

5. 5. 6. 2 Data  Flow 

• Input  Data 

RR(ifj)  - range-ordered  list  of  marks  (array  RRDOT) 

• Output  Data 

RR(i,j)  - range-ordered  list  of  marks  (without  redundant 
marks)  (array  RRDOT) 

5. 5. 6. 3 Detailed  Description 

REMOVE  examines  each  type  1 mark  (true  peaks),  searchs  for  and 
deletes  any  adjacent  type  2 marks.  Since  the  list  is  range-ordered, 
each  search  should  involve  only  a few  adjacent  marks.  Deleted  marks 
are  purged  from  the  list. 

The  detailed  flowchart  of  REMOVE  is  shown  in  Figure  5.8. 
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6.0  COINCIDENCE  DETECTION  ALGORITHM 


6.1  PURPOSE  OF  SUBSYSTEM 

The  primary  purpose  of  coincidence  detection  is  the  removal  of  marks 
caused  by  range  ambiguities  of  the  burst  waveform. 

Coincidence  detection  processing  is  required  only  for  those  trans- 
missions on  which  a burst  pair  is  employed. 

6 . 2 GENERAL  DESCRIPTION 

The  coincidence  detection  algorithm  operates  on  the  range,  range-rate 
marks  and  the  sampled  video  from  two  bursts  on  the  same  transmission  to 
produce  a set  of  marks  for  that  transmission.  In  the  bulk  filtering 
application  described  here,  a novel  implementation  of  coincidence  detection 
is  used  which  has  two  equally  important  objectives 

• Removing  marks  caused  by  range  ambiguities 
of  the  burst  waveform 

e Achieving  a high  probability  of  RV  detection  in  the 
presence  of  interference  by  fragments  and  fragment 
ambiguities . 

This  design  has  the  advantage  that  in  an  environment  with  severe  interfer- 
ence the  algorithm  sacrifices  ambiguity  removal  performance  to  maintain 
an  acceptable  probability  of  detection,  while  in  less  severe  interference 
the  ambiguity  removal  performance  improves  substantially. 

The  coincidence  detection  algorithm  recommended  for  the  coherent 
bulk  filter  is  designed  to  give  equal  Importance  to  the  removal  of  marks 
caused  by  range  ambiguities  and  the  minimization  of  RV  leakage.  To 
accomplish  this,  the  algorithm  departs  from  the  conventional  method  of 
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Implementing  coincidence  detection  which  requires  that  marks  from  the 
two  bursts  occur  within  a specified  distance  in  range  and  range-rate. 

The  recommended  algorithm  is  defined  as  follows: 

A coincidence  detection  occurs  at  a given  range  and 
range-rate  if  (1)  a range,  range-rate  mark  is  generated 
on  either  of  the  two  bursts  and,  (2)  the  video  from  the 
other  burst  exceeds  a threshold. 

This  algorithm  reduces  RV  leakage  in  the  important  situation  where  a 
large  fragment  ambiguity  is  aligned  in  range  with  the  RV  return  on  one 
of  the  bursts.  Since  the  fragment  range-rate  is  lower  chan  chat  of  the 
RV,  a mark  is  generated  only  at  the  lower  range-rate.  On  the  other  burst, 
Che  ambiguity  is  shifted  by  virtue  of  the  different  pulse  spacing,  and  a 
mark  is  generated  at  the  RV  range-rate.  The  recommended  algorithm  then 
generates  a coincidence  detection  from  Che  mark  on  the  burst  without 
Interference  and  the  video  above  the  threshold  on  Che  burst  without  inter- 
ference. The  conventional  implementation  of  coincidence  detection  would 
not  generate  a detection  for  this  case. 

An  overview  flowchart  for  the  coincidence  detection  subsystem  is 
shown  in  Figure  6.1.  As  shown  in  the  overview  flowchart,  the  coincidence 
detection  algorithm  first  compares  the  marks  from  burst  number  1 with  the 
video  from  burst  number  2,  and  Chen  compares  the  marks  from  burst  number  2 
with  the  video  from  burst  1.  The  two  lists  of  marks  are  then  merged  and 
sorted.  Finally,  redundant  coincidence  detections  on  Che  same  object, 
resulting  from  marks  on  both  bursts,  are  Identified  and  removed. 

The  input  data  consists  of  marks  from  the  marking  algorithm  and 
video  data.  The  output  consists  of  a list  of  marks  passing  coincidence 
detection.  This  list  is  passed  to  the  designation  sequence  combining 
algorithm. 
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Figure  6.1  Overview  Flowchart  for  Coincidence  Detection 
Subsystem 
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6.3  DATA  SPECIFICATION 

1 

i 

6.3.1  Input  Data 

1 

Program 

Symbol 

Engineering 

Symbol 

Description  I 

RRDOTi(Q,.M) 

Ri^(q,m) 

1 

Range-ordered  list,  indexed  by  q,  of 
range,  range-rate  marks  from  burst  I 

number  1 specified  by  range  (m=l,m) , • 

range-rate  (m-2,m/s),  and  a flag  (m=3) 
used  during  marking.  j 

NMl 

NM^ 

Number  of  range,  range-rate  marks  for 
burst  number  1. 

RRDOT2(Q,M) 

RR2(q,m) 

Range-ordered  list,  indexed  by  q,  of 
range  range-rate  marks  from  burst 
number  2 specified  by  range  (m=l,m), 
range-rate  (m-2,m/s),  and  a flag  (m*3) 
used  during  marking. 

NM2 

NM2 

Number  of  range,  range-rate  marks  for 
burst  number  2. 

t 

RS 

^S 

Range  corresponding  to  start  of  video 
data  processing  (m) . 

DELR 

A 

Range  sample  spacing  (m) . 

DS 

°S 

Range-rate  corresponding  to  lowest 
Doppler  channel  for  which  video  data  is 
processed  (m/s) . 

DELD 

Range-rate  sample  spacing  (Doppler  channel 
spacing)  (m/s). 

VIDEOKK.I) 

rj^Ot.i) 

Amplitude  of  video  at  range  R +(k-l)A 
and  range-rate  D^+-(i+l)C  for  Burst 
number  2 (ra) . 

RREF 

^0 

Reference  range  for  specification  of 
video  threshold  (m) . 

THRCD 

^CD 

Video  threshold  (specified  as  a radar 
cross  section)  at  the  reference  range 
R used  to  generate  a coincidence 
detection  (m2). 
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6.3.1  (Continued) 


Program 

Svmbol 


RDOTCD 


Engineering 

Svmbol 


Description 

Rai.ge  interval  used  in  the  removal 
of  redundant  marks  on  the  same  object 
at  the  end  of  the  coincidence  detection 
operation  (ra) . 

Range-rate  interval  used  in  the  removal 
of  redundant  marks  on  the  same  object 
at  the  end  of  the  coincidence  detection 
operation  (m/s). 


6.3.2  Output  Data 


CDMRKS(L,M)  RRR(S.,m) 


Range-ordered  list,  indexed  by  i, 
of  range,  range-rate  marks  that  have 
passed  coincidence  detection. 

m =■  1 is  range  (m) 
m ” 2 is  range-rate  (m/sec) 
m = 3 is  angular  coordinate,  u 
m » 4 is  angular  coordinate,  v 
m = 5 is  the  burst  number. 

Number  ot  range,  range-rate  markts 
that  have  passed  coincidence  detection. 


6.3.3  Parameter  Settings 


Four  parameters  can  be  specified  in  the  coincidence  detection  algorithm. 
The  reference  range,  R^,  at  which  the  video  threshold  is  specified  is  normally 
the  middle  of  the  range  extent  corresponding  to  the  acquisition  band,  i.e., 

Rq  ■ 1/2  (Rgy  video  threshold,  which  is  applied  to  the 

burst  opposite  that  containing  the  mark,  is  set  to  achieve  the  desired  operat- 
ing point  (Pjj,  at  the  output  of  coincidence  detection.  Only  cases  with 

the  video  threshold,  set  equal  to  the  detection  threshold,  T^,  have 

been  simulated.  This  was  found  to  be  a suitable  operating  point,  and 
is  recommended.  The  other  two  parameters,  R^^^  and  the  range  and  range- 

rate  intervals,  respectively,  used  in  the  removal  of  redundant  marks  on  the 
same  object,  are  selected  based  on  the  resolution  properties  of  the  waveform. 
It  is  recommended  that  they  be  set  equal  to  twice  the  measurement  covariances, 
i.e.  R^jj  - 2Cj^  and  R^  - 2a^. 


6.4  MATHEMATICAL  RELATIONSHIPS 


A general  flow  diagram  for  the  coincidence  detection  algorithm  is 
shown  in  Figure  6.1. 


The  first  step  in  the  processing  compares  the  range,  range-rate  marks 
from  burst  number  1 with  the  video  from  burst  number  2.  For  each  mark, 
RR^(i,  m) , from  burst  number  1,  the  indices  of  the  closest  video  sample  of 
burst  number  2 are  determined.  The  index  in  range  is 


k - Integer 


(6.1) 


and  the  index  in  range-rate  is 


1 ■ Integer 


- Dg 


(6.2) 


The  video  amplitude  of  burst  number  2 at  range  sample  k and  range-rate 
sample  i is  then  compared  with  the  video  threshold,  TCj^,  given  by 


If 


0 

k - 1)  A 


2 


r.(k,l)  > TCj^  . 


(6.3) 


the  range  range-rate  mark  RR^(i,m)  is  retained;  otherwise,  it  is  deleted. 
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Similarly,  the  range,  range-rate  marks  from  burst  number  2 are 
compared  with  the  video  from  burst  number  1.  For  each  mark,  m) , 

the  closest  video  sample  from  burst  number  1 has  range  index 


k - meager  ~ ^ il 

I 1 2 ) 


and  range-rate  index 

1 “ Integer  ^^'1 


RR, (1,2)  - D_  3 ) 
i 


(6.4) 


(6.5) 


If 


r^(k,i)>TCj^  , 


the  mark  ^^(l,  m)  is  retained;  otherwise,  it  is  deleted. 

The  range,  range-rate  marks  retained  from  both  bursts  are  merged 
into  a range-ordered  list,  RRR  (q,m),  indexed  on  q,  of  marks  specified 
by  range  (m  - 1) , range-rate  (m  - 2),  and  burst  number  (m  - 3)  on  which 
the  mark  was  generated.  A test  is  then  performed  to  identify  redundant 
marks  on  the  same  object  caused  by  the  occurrence  of  marks  on  both  bursts, 
If  two  marks  RRR(Z,m)  and  RRR(q,m)  are  found  for  which 


RRR(i,3)  ^ RRR(q,3) 


(6.6) 


and 


|RRR(i,l)  - RRR(q,l)  \ 1 R, 


CD 


(6.7) 


and 


RRR(i,2)  - RRR(q,2)  £ R^ 


(6.8) 


I 


the  mark  with  the  minimum  range  is  retained,  and  the  other  mairk  la 
deleted.  (This  is  an  arbitrary  choice  of  which  mark  to  retain.) 
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6.5  DETAILED  DESCRIPTIONS 


I 

I 

6.5.1  Executive  Subroutine  | 

6. 5. 1.1  Purpose  | 

The  coincidence  detection  executive  subroutine  controls  the  data  | 

flow  and  the  execution  of  subroutines  for  this  subsystem. 

6. 5. 1.2  Detailed  Description 

The  operation  of  the  executive  subroutine  is  straightforward.  Sub- 
routins  are  called  as  shown  in  the  overview  flowchart  in  Figure  6.1. 


i 

' 
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6.5.2  Subroutine  COIN 

6. 5. 2.1  Purpose 

For  each  mark  from  a burst,  COIN  checks  for  a threshold  crossing  of 
the  video  from  the  other  burst.  If  the  video  does  not  exceed  the  threshold, 
the  mark  is  rejected. 

6. 5. 2. 2 Essential  Input/Output  Data 


Range,  range-rate  marks  from  this  burst 
(burst  n) 

Video  from  the  other  burst  (burst  n') 


Range,  range-rate  marks  from  this  burst 
pass  first  step  of  coincidence  detection 

6. 5. 2. 3  Detailed  Description 

Subroutine  COIN  examines  each  mark  from  a burst.  The  range,  range- 
rate  sample  closest  to  the  mark  is  determined.  The  video  of  the  other 
burst  is  examined  at  that  sample  to  check  for  a threshold  crossing.  A 
detailed  flowchart  for  COIN  is  shown  in  Figure  6.2. 


I 

I 


6.5.3  Subroutine  MERGE 


6. 5. 3.1  Purpose 

>IERGE  combines  the  marks  from  the  first  and  second  bursts  into 
one  list,  and  then  sorts  the  list  in  range  order.  Marks  which  have  been 
flagged  for  deletion  may  be  rejected  at  this  point.  The  merge  and  sort 
operations  are  self  explanatory. 
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6.5.4  Subroutine  CULL 


6. 5. 4.1  Purpose 

CULL  removes  redundant  marks  caused  by  the  same  object. 

6. 5.4.2  Essential  Input/Output  Data 
Input 


RRDOT(q,m)  - range-ordered  list  of  marks  from 
both  bursts 

Output 

RRDOT(q,m)  - range-ordered  list  of  marks  which  passed 
test  with  redundar.t  marks  removed 

6. 5.4. 3  Detailed  Description 

For  each  mark,  CULL  examines  adjacent  marks  over  a range  extent  of 

+R  ^ . Only  marks  from  the  other  burst  are  examined.  If  another  mark 
-CD  ^ 

is  found  whose  range-rate  differs  from  the  current  mark's  by  less  than 
^CD  ’ ®ark  with  the  lesser  range  is  deleted. 


A detailed  flowchart  for  CULL  is  shown  in  Figure  6.3. 


Figure  6.3  Subroutine  CULL  Detailed  Flowchart 


I 
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7.0  COHERENT  DOUBLE  GATING  (CDG)  ALGORITHM 


7.1  PURPOSE  OF  SUBSYSTEM 

The  coherent  double  gating  algorithm  performs  two  general  functions: 
(1)  the  removal  of  detections  caused  by  clutter  and  ambiguities  chat  are 
not  removed  by  Che  single  transmission  processing,  and  (2)  the  generation 
of  a state  estimate  on  each  threatening  target  for  track  initiation.  The 
CDG  algorithm  performs  the  same  processing  after  each  transmission, 

it 

regardless  of  whether  it  is  a burst  pair  or  a single  burst. 

7.2  GENERAL  DESCRIPTION 

The  CDG  algorithm  consists  of  three  tests  for  reducing  the  number  of 
potentially  threatening  objects  designated  in  the  beam  and  a maximum 
likelihood  estimation  computation  for  generating  the  best  state  estimate 
for  track  initiation. 

The  first  test  is  a deceleration  or  gating  test  which  eliminates  any 
object  whose  trajectory  deviates  by  some  amount  from  a constant  range-rate 
trajectory.  To  perform  this  test,  the  maximum  likelihood  estimate  of  the 
two-element  state  (R,R)  is  computed  for  each  object  after  transmission  j-1. 
Based  on  this  estimate  and  assuming  a constant  range-rate  trajectory,  an 
uncertainty  region  (elliptical  gate)  is  computed  in  the  range,  range-rate 
plane  for  the  measurement  on  transmission  j . If  the  measurement  (mark) 
falls  within  the  gate,  the  object  is  retained.  If  the  measurement  falls 
outside  of  the  gate,  a miss  is  recorded.  If  the  number  of  allowable  misses 
is  exceeded,  the  object  is  dropped  from  further  consideration.  If  the 
number  of  allowable  misses  is  not  exceeded,  the  object's  state  estimate  is 
updated  to  the  present  transmission  time.  Since  a miss  can  be  caused  by 
failure  to  detect  the  object,  it  is  preferable  to  allow  at  least  one  miss 

Although  the  CDG  algorithut  has  tuis  CApa'oility , this  report 
specifies  Its  use  only  with  burst  pairs. 
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during  Che  sequence  of  transmissions  co  achieve  a high  probability  of  desig- 
nation at  the  end  of  the  sequence  without  requiring  an  unreasonably  high 
single-hit  probability  of  detection.  A strategy  allowing  M misses  in 
a total  of  N transmissions  is  denoted  by  N-M:N. 

The  CDG  algorithm  uses  one  of  two  methods  to  achieve  an  N-M:N  strategy 
which  yield  different  operating  points.  The  type  1 N-M:N  strategy,  or 
N-M:N(1),  requires  chat  the  object  be  detected  on  the  first  transmission, 
i.e.,  no  new  object  state  estimates  will  be  initiated  on  transmissions  2 
through  N.  The  type  2 N-M:N  strategy,  or  N-M:N(2),  allows  a miss  to  occur 
on  the  first  transmission.  This  means  chat  marks  not  falling  into  an  ellip- 
tical gate  on  the  second  transmission  are  used  co  initialize  new  object  state 
estimates .with  one  miss.  The  type  2 strategy  allows  entries  to  be  added 
to  the  CDG  file  on  two  transmissions,  and  hence  achieves  a higher  proba- 
bility of  designation  at  the  expense  of  more  computations  and  a slightly 
higher  level  of  ghosting. 

The  second  test  in  the  CDG  algorithm  is  a velocity  test.  This  test 
utilizes  the  increased  confidence  in  the  range-rate  estimate  from  multiple 
transmissions  to  effectively  sharpen  the  cut-off  at  the  edges  of  the  RV 
range-rate  acceptance  region.  It  is  used  after  each  transmission  beyond 
Che  first  and  tests  whether  the  range-race  estimate  is,  with  high  confidence, 
outside  of  the  acceptance  region. 

The  third  test,  called  an  association  test,  is  designed  to  eliminate 
multiple  entries  in  the  CDG  algorithm  file  corresponding  to  the  same  object. 
Such  multiple  entries  can  arise  from  multiple  marks  on  the  same  object  from 
the  first  transmission,  or  from  multiple  marks  falling  in  the  gate  for  Che 
gating  test  on  the  second  transmission.  The  association  test  is  performed 
after  each  transmission  beyond  the  second. 

The  maximum  likelihood  estimation  of  (R,R)  for  crack  initiation  is 
generated  as  part  of  the  deceleration  or  gating  test. 
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Figure  7.1  Is  an  overview  flowchart  for  the  coherent  double  gating 
f algorithm.  During  the  first  iteration  shown  in  the  flowchart  (j"l)  the 

list  of  object  state  estimates  is  initialized  with  the  marks  from  trans- 
mission number  one.  These  marks  can  be  generated  by  the  processing 
described  either  in  Section  6 for  burst  pair  transmissions  or  in  Section  5 
for  single  burst  transmissions. 

On  the  second  iteration  ( j»2) , each  object  state  estimate  from  iteration 
one  is  used  to  establish  a range,  range-rate  gate  for  the  deceleration 
test.  Marks  from  transmission  two  passing  the  gating  test  are  used  to 
update  state  estimates  of  the  object.  Those  objects  which  then  pass  the 
velocity  test  are  placed  in  the  updated  object  file.  Finally,  if  the 
strategy  type  is  2,  marks  not  falling  in  any  gate  are  used  to  initiate 
new  object  state  estimates. 

From  the  third  iteration  on,  state  estimates  are  updated  by  the  mark 
closest  to  the  center  of  the  gate.  Each  object  must  then  pass  the  velo- 
city test.  The  remaining  objects  are  processed  by  the  association  test 
to  -emove  redundant  tracks.  After  the  last  iteration,  the  list  of  objects 
is  passed  on  to  the  Known  Object  Recognition  Algorithm. 


FOR  EACHJRANSMISSiON, J_ 


COGALG 


j-1? 


MARKS  FROM  I 
TRANSMISSION 


PREVIOUS  STATE 

msi^ 


INIT  ^ 

initialize  (56Jcc1 

STATE  ESTIMATES/ 


UF'^ATE 2 

GAiING  TEST.  STATE  UPDATE , 
VELOCITY  TEST,  AND  STRATEGY/ 
TYPE  IMPLEMENTATION  / 


1 CURPENT  STATE | 
\ESTIMATE  OF  OBJECTS 


> \ 

\ \ 

\ \ NO 

\ \ 

\ \ 

\ \ 


j >_  3? 


RETURN 


REMOVE  REDUNDANT 
STATE  ESTIMATES 


Figure  7.1  Overview  Flowchart  for  Coherent  Double 
Gating  Algorithm  Subsystem 
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7.3  DATA  SPECIFICATION 


7.3.1  Inrut  Data 


The  following  data  is  input  via  the  CDG  OPERATING  PARAMETERS  File: 


Program  Engineering 

Symbol  Symbol Descripcion 


NTRAN 

NMISS 

IS 


CG 


CV 

RDTMIN 

RDTMAX 

CA 


N 


Number  of  transmissions  in  CDG  sequence 


M 

IS 


CG 


CV 


min 


R 

max 

CA 


Number  of  misses  allowed  La  CDG  sequence 

Specifies  an  N-M:N(IS)  algorithm  (IS=*1 
requires  a detection  on  the  first  trans- 
mission, IS»2  requires  a detection  on 
first  or  second  transmission) . 

Constant  related  to  the  probability  of  a 
detection  from  a threatening  target  falling 
In  the  elliptical  gate. 

Constant  related  to  the  probability  of  a 
threatening  target  passing  the  velocity  test. 

Lower  range-rate  limit  in  velocity  test  (m/sec) 

Upper  range-rate  limit  in  velocity  test  (m/sec) 

Constant  related  to  the  probability  of 
incorrectly  associating  two  detections  from 
different  objects. 


The  following  data  is  input  via  the  COINCIDENCE  DETECTION  MARKS  File: 


CDMRKS(L,M)  RR  (i,m)  Range-ordered  list,  indexed  by  q,  of  range, 

range-race  marks  from  transmission  j 

m ■ 1 range  (m) 
m-2  range-rate  (m/sec) 
m "3  angular  coordinate,  U 
m “4  angular  coordinate,  V 
tn  "5  burst  number 
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NCD  NM^ 

IMEAS  t . 

J 


Number  of  range,  range-rate  marks 
in  array  CDMRKS 

Time  of  transmission  j (sec) . 


The  following  data  is  input  via  the  CDG  TRACK  file: 


NTN 


CDGTRK(K.M) 


NT 


j-1 


Number  of  objects  passing  the  CDG  Algorithm 
on  transmission  j-1  (number  of  entries 
1»  DT^_p  . 


DT._j^(k,m)  Range-ordered  designation  list,  indexed 
by  k,  of  objects  remaining  in  the  CDG 
algorithm  after  processing  of  transmission  j. 


M = 1 Range  estimate  (m) 

M = 2 Range-rate  estimate  (m/s) 

M = 3 Number  of  misses 

2 

M = 4 Variance  of  range  estimate  (m  ) 

M = 5 Covariance  of  range  and  range -rate 

estimates  (m'^/s) 

M =>  6 Variance  of  range-rate  estimate  (m‘’/s"‘) 

M = 7 Angular  coord:' nate  U 

M = 8 Angular  coordinate  V 


The  following  data  is  input  via  the  MEASUREMENT  PARAMETERS  file: 


SIGR  o 
SIGRDT  a 


A-priori  range  error  standard  deviation  (m) 

A-priori  range-rate  error  standard  deviation 
(m/ sec) 


7.3.2  Ouput  Data 

The  following  data  is  output  via  the  CDG  TRACK  file: 


NTN 


CDGTRK  (K,M) 


NT.  Number  of  objects  oassing  the  CDG  algorithm 

•*  on  transmission  j 

DT,(k,m)  Updated  designation  list  of  objects 

^ remaining  after  processing  transmissio.-  j 
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7.3.3  Parameter  Settings 


The  ten  parameters  N,  M,  IS,  a , a' , CG,  CV,  R ^ , R , and  CA  are 

r r^  _ min  max 

used  exclusively  in  the  03G  Algorithm.  R . and  R are  set  to  the 

min  max 

minimum  and  maximum  range-rates  of  RV's  threatening  the  defended  area. 
Values  of  -7.2  km/s  and  -5.4  km/ s respectively  were  used  in  all 
simulations.  The  constant  CV  is  set  to  2.57  to  achieve  a probability  of 
.995  of  a threatening  target  passing  the  velocity  test.  Simulation 
reSL.ts  with  bursts  of  16  10  MHz  bandwidth  pulses  and  coincidence  detec- 
tion indicate  tlu.  t settings  of 

N-M:N(IS)  - 5:6(2),  4:6(2),  or  4:5(2) 

(a  ,0")  * 3.66  m,  156  m/s 

r r 

CG  =»  10.6 

CA  » 5.0  to  10.0 

are  suitable  operating  points.  For  32  pulse  bursts  and  can  be 
reduced  to  approximately  1.83  m and  38.9  m/s. 


7.4  MATHEMATICAL  RELATIONS 


A general  flow  diagram  for  the  CDG  algorithm  is  shown  in  Figure  7.1. 


After  the  first  (search)  transmission,  j=l,  variables  are  initialized 


as  follows. 


ND^  - NMj 


(7.1) 


For  each  q » l,2,...,NMj,  set 


DT^(q.l)  - RRj(q,l) 
DT^(q,2)  = RRj(q,2) 
DTj(q,3)  = 0 
DT^(q,4)  - 
DT^(q,5)  = 0 
DT^  (q,6)  =•  0^ 


(7.2) 


The  remaining  operations  are  performed  after  the  second  (first  verify) 
and  subsequent  transmissions,  j*2,3,...  First,  the  measurement  error 
covariance  matrix  is  set  to 


(7.3) 


the  transition  matrix  is  set  to 


1 


1 

^ J 


(7.4) 


and  the  count  of  objects  tn  DT 


j 


is  set  to  zero 


NDj  - 0 . 


(7.5) 


For  each  object  tn  the  CDG  algorithm  after  the  processing  of  transmission 
j-1,  that  is  for  each  k - 1,2, , . . the  following  sequence  of  computa' 

tions  is  performed. 


The  state  estimate  for  object  k after  transmission  j-1  is 
set  to  _ 


DT^_^(k,l) 

DTj_^(k,2) 


and  the  estimate  covariance  matrix  for  object  k is  set  to 


DTj_i(k,4)  DT^_^(k,5) 


L 


DT^_^(k,5)  DTj_j^(k,6) 


The  predicted  state  for  object  k at  time  is 


(7.6) 


and  the  predicted  covariance  for  object  k at  time  t^  is 


P,  • <P  P,  . 
k k 


(7.7) 


The  covariance  matrix  for  the  difference  between  the  predicted  state 
for  object  k and  any  one  of  the  marks  from  transmission  j is 


\ - ^k  " 


(7.8) 
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For  each  mark  q,  q*l, 2, . , . ,NM^ , from  transmission  j,  the  difference 
between  the  predicted  state  for  object  k and  the  measured  state  object 


q is 


I RRj(q,l) 

- I 


RRj(q,2) 


(7.9) 


The  mark  q falls  within  the  acceptance  gate  of  object  k if 


At  this  point,  two  different  paths  are  identified:  Path  I,  if  no  narks 
fall  within  the  acceptance  gate;  and  Path  II,  if  one  or  more  narks  fall 
within  the  acceptance  gate. 


Path  I:  No  Marks  in  Acceptance  Gate 


If 

DTj_^(k,3)  - M , (7.11) 

where  M is  the  allowed  number  of  misses,  delete  object  k from  further 
tracking  by  simply  beginning  processing  on  object  k+1.  Object  k will  not 
appear  In  the  updated  BEAM  DESIGNATION  FILE  in  that  case. 

If 


DTj_j^(k,3)  < M , 

increment  the  designation  counter 

NDj  - ND^  + 1, 

increment  tnc  ’’auss"  counter  for  ooject  k 

DT,(ND  ,3)  - DT,  ,(k,3)  + 1 
J j J"-*- 


(7.12) 


(7.13) 


7-10 


and  set  the  other  elements  of  the  BEAM  DESIGNATION  FILE  DTj (ND^ ,m) , equal 
to  the  predicted  state  and  covariance  of  object  k: 


DT^(ND^,1)  - 

DTj(ND^,2)  - 

X^(2) 

DT^(ND^,4)  * 

P^(l,l) 

DTj(ND^,5)  - 

P,(1.2) 

DTj(NDj,6)  - 

Pj^(2,2) 

Path  II:  Marks  in  Acceptance  Gate 

Within  Path  II,  three  cases  can  arise: 

(1)  If  j > 2 and  more  chan  one  detection  falls  within  the 

acceptance  gate  of  object  k,  only  the  mark  q with  the 

T -1 

smallest  value  of  Z,  V,  Z,  is  retained  and  the  others 
Tt,q  k k,q 

are  deleted. 

(2)  If  j * .2  and  more  chan  one  mark  falls  within  the 
acceptance  gate,  all  such  marks  are  retained  for 
further  processing.  If  j*2,  the  index,  q of  each 
mark  retained  is  stored  in  the  array,  J.  This  array 
is  required  to  implement  an  N-M:N(2)  strategy. 

(3)  If  one  mark  falls  within  the  acceptance  gate,  it 
is  retained  for  further  processing. 


For  each  mark  retained  by  the  gating  test,  the  state  estimate  which 
results  from  the  combination  of  object  k and  mark  q 


X 

q 


(7.15) 
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where 


S.  ■ Vk' 

and  Che  estimate  covariance  matrix 

- (I  - K^)P^  (7.16) 

are  updated,  where  I is  the  identity  matrix.  Each  of  these  objects  is 
processed  through  the  velocity  test: 


R + CV  * \/p(2,2)  > X (2)  > R , ~ (y  * n/p„(2,2)  (7.17) 

max  q — q ~ min  q . 

If  the  object  passes  this  test,  then  the  object  is  retained  in  the 

updated  BEAM  DESIGNATION  FILE;  otherwise,  Che  object  is  deleted.  For 

each  object  retained,  set 


and 


NDj  - NDj  + 1 


DTj(NDj,l) 

DTj(NDj,2) 

DTj(NDj,3) 

DTj(NDj,4) 

DTj(ND^,5) 

DTj(NDj,6) 


- 

- X (2) 

q 

- DTj_^(k,3) 

- 

- Pjl.2) 

i 

- PJ2,2) 

q 


(7.18) 


(7.19) 


This  is  the  end  of  the  processing  for  each  object,  k,  which  remained 
in  Che  CDG  algorithm  after  transmission  j-1. 
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If  a type  2 strategy  is  used,  the  extra  steps  described  in 
this  paragraph  must  be  carried  out  for  J*2,  i.e.,  after  the  second 
transmission.  Each  mark,  not  falling  within  an  acceptance  gate,  i.e., 
q J,  is  stored  in  the  list  of  objects  retained  by  the  CDG  Algorithm 
after  transmission  j.  Each  is  treated  as  having  been  "missed" 
on  the  first  transmission.  For  each  such  mark,  q,  set 


ND^  - NDj  + 1 (7.20) 

and 

DT^(NDj,l)  - RRj(q,l) 

DT^(NDj,2)  « RR^(q,2) 

DTXND^,3)  - 1 

DT,(ND, ,4)  - a2 
J J R 

DTj(NDj,5)  - 0 

DT^(ND, ,6)  - a?  (7.21) 

i i R 


When  the  list  DTj(l,m)  is  completed,  all  entries  are  range-ordered. 


For  each  transmission  beyond  the  second,  j >2,  an  association  test 
is  performed  on  all  objects  retained  in  the  CDG  algorithm  to  eliminate 
redundant  entries  from  the  same  object:  two  entries,  i and  m,  are  said 

to  be  from  the  same  object  (associated)  if 

D 1 CA  (7.22) 

where 


DTj(l,4) 

DTjU,5) 

DTj(m,4) 

DT^(m,5) 

B - 

DT  <'J.,5) 

L ^ ' 

DT 

+ 

DTj(m,5) 

DTj(m,6) 
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If  entries  I and  m are  associated,  the  entry  with  the  greater  number 
of  misses  is  deleted  or,  if  DTj(i,3)  and  DT^(m,3)  are  equal,  the  entry 
with  the  greater  range  is  deleted. 


7 . 5 DETAILED  DESCRIPTIONS 


7.5.1  Executive  Subroutine 

i 

7. 5. 1.1  Purpose 

The  coherent  double  gating  algorithm  executive  subroutine  controls 
I data  flow  and  execution  of  subroutines  for  this  subsystem. 

7. 5. 1.2  Detailed  Description 

The  executive  subroutine  calls  subroutines  as  shown  in  the  overview 
^ flowchart.  Figure  7.1. 

r 

f 

I 

I 

I 

I 

I 

I 

I 
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7.5.2  Subroutine  INIT 


7. 5. 2.1  Purpose 

On  the  first  transmission,  INIT  performs  the  initialization 
described  by  equations  (7.1)  and  (7.2).  The  array  DT^^  is  then  ready 
for  update  on  subsequent  transmissions.  A detailed  flowchart  of 
subroutine  INIT  is  shown  in  Figure  7.2. 
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7.5.3  Subroutine  UPDATE 


7. 5. 3.1  Purpose 

Subroutine  UPDATE  propagates  state  estimates  from  the  previous 
transmission  time  to  the  present  transmission  time  and  implements  the 
strategy  type.  UPDATE  calls  subroutine  TESTS  to  perform  the  velocity 
test  on  updated  state  estimates. 

7. 5. 3. 2 Essential  Input/Output  Data 
Input 


RRj(q,i)  range,  range-rate  marks  from  j transmission 
DTj  l(k,m)  object  file  from  previous  transmission  (j-1) 


Output 


DTj(k,m)  current  object  file 


7.5. 3.3  Detailed  Description 


For  each  object  in  DT^  j^(k,m)  , UPDATE  performs  the  processing 
described  by  equations  (7.3)  through  (7.14)  and  by  equations  (7.20) 
through  (7.21).  A detailed  flowchart  of  UPDATE  is  shown  in  Figure  7.3. 
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I 

I 


Figure  7.3  Subroutine  UPDATE  Detailed  Flowchart 
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7.5.4  Subroutine  TESTS 


7. 5.4.1  Purpose 

Subroutine  TESTS  performs  two  types  of  tests  on  the  marks  which  fall 
within  the  acceptance  region  of  an  object’s  updated  state  estimate.  These 
tests  involve  the  choice  of  objects  to  be  retained,  and  the  velocity  test. 


7. 5.4. 2 Input/Output  Data 


Input 


Program 

Symbol 

ZQK(I) 

K 

VKINV(2,2) 

XK(2) 

PK(2,2) 

RDTMIN 

RDIMAX 

CV 


Engineering 
Symbol 


^ittin 

^max 

CV 


Description 

List  of  the  indices,  q,  of  each  mark 
in  REDOT  array  which  falls  in  acceptance 
region  of  object  k 

List  of  differences  between  predicted 
and  measured  values  for  each  mark  in 
array  QMARKS(I) 

Index  in  DT^^j^Ck.m)  object  file  for 
current  object 

Inverse  error  covariance  matrix  for  object  k 

Predicted  state  for  object  k 

Predicted  covariance  for  object  k 

Lower  range~rate  limit  for  velocity  test 

Uppex  range~rate  limit  for  velocity  test 

Constant  related  to  the  probability  of 
a threatening  target  passing  the 
velocity  test 


7. 5. 4. 2 (Concinued) 


Figure  7.4  Subroutine  TESTS  Detailed  Flowchart 
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7.5.5  Subroutine  ASSOC 


7. 5. 5.1  Purpose 

The  purpose  of  ASSOC  is  the  removal  of  multiple  entries  in  the  BEAM 
DESIGNATION  FILE  arising  from  the  same  object. 

7. 5. 5. 2 Input /Output  Data 


Input 

Program 

Engineering 

Symbol 

Symbol 

DESVEC(I.K) 

DTj(m,k) 

Output 

DESVEC(1,K) 

DTj(m.k) 

7. 5. 5. 3  Detailed  Description 


Description 

List  of  objects  after  transmission  j. 

List  of  objects  after  transmision  j with 
redundant  tracks  removed. 


Each  pair  of  objects  retained  by  the  CDG  Algorithm  is  tested  to 
determine  whether  the  objects  are  associated.  If  they  are  associated, 
then  the  object  with  the  greater  number  of  misses  is  deleted  unless  both 
objects  have  an  equal  number  of  misses  in  which  case  Che  one  with  Che 
greater  range  is  deleted. 

ASSOC  performs  Che  processing  described  by  equations  (7.22)  through 
(7.24).  The  detailed  flowchart  of  ASSOC  is  given  in  Figure  7.5. 
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Figure  7.5  Subroutine  ASSOC  Detailed  Flowchart 


8.0  OBJECT  BEAM  POSITION  ALGORITHM 


8.1  PURPOSE  OF  SUBSYST01 

The  object  beam  position  algorithm  computes  the  UV  coordinates 
for  designated  objects  and  stores  those  coordinates  In  the  BEAM 
DESICaiATION  FILE. 

8.2  GENERAL  DESCRIPTION 

Two  channels  of  monopulse  data  are  available  which  provide  object 
angular  position  (AU  and  AV)  with  respect  to  the  beam  center. 

The  object  beam  position  algorithm  obtains  object  angular  measure- 
ments from  the  monopulse  data  for  each  object  that  passes  coincidence 
detection.  These  monopulse  measurements  are  combined  with  the  beam  UV 
coordinates  to  obtain  object  UV  coordinates. 

The  overview  flowchart  for  the  Object  Beam  Position  Algorithm  is 
shown  In  Figure  8.1. 
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COMPUTE 

OBJECT  ANGULAR 
POSITION 


Figure  8.1  Object  Beam  Position 


Algorithm  Overview  Flowchart 


8.3  DATA  FLOW 


8.3.1  Input  Data 

The  following  data  is  input  from  the  ALF  via  the  Beam  Status  file: 


Program 

Symbol 

ID 

IV 


Engineering 
Symbol 


Description 
Beam  U index 
Beam  V index 


The  following  data  is  input  from  the  Coincidence  Detection  Algorithm 
via  the  COINCIDENCE  DETECTION  MARKS  file. 


MCD  - Number  of  coincidence  detection  marks 

from  the  current  burst  pair 

CDMRKS(L,M)  RRR(i,m)  Range,  range-rate  marks  from  coincidence 

detection 


R 

m ■ 1 

is 

range  (m) 

R 

m ” 2 

is 

range-rate  (m/sec) 

K 

m • 5 

is 

the  burst  nximber 

I is  the  mark  index 


The  following  data  is  input  from  the  Matched  Filter  Processing  Function 
via  the  BURST  1 MONOPULSE  file  and  the  BURST  2 M(WOPULSE  file: 


RSM 

REM 

DEL 

DSM 


*4 


Range  corresponding  to  start  of  video 
data  processing  (m) 


Range  corresponding  to  end  of  video 
data  processing  (m) 


L 


D 


S 


P 
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Range  sample  spacing  (m) 

Range-rate  corresponding  to  lowest 
Doppler  channel  for  which  video 
data  is  processed  (m/s) 


-rtysiiTjir- 


8.3.1  (Continued) 


DEM 

“e 

Range-rate  corresponding  to  highest 
Doppler  channel  for  which  video  data 
is  processed  (m/s) 

DEED 

Range-rate  sample  spacing  (Doppler 
channel  spacing)  (m/s) 

UBST1(K,L) 

“bj 

Monopulse  U measurement,  beam  1 

VBST1(K,L) 

Monopulse  V measurement,  beam  1 

UBST2(K,L) 

< 

Monopulse  U measurement,  beam  2 

VBST2(K,L) 

Monopulse  V measurement,  beam  2 

For  the  four  monopulse  arrays,  note  that  K - 1 at  and  L ■ 1 at  . 

The  following  data  is  input  via  the  OBJECT  BEAM  POSITION  PARAMETERS 

file: 


UBEAM(IO,IV)  0 U coordinate  (sine  space)  of  the 

° lU,  IV  beam 

VBEAM(IU,IV)  V V coordinate  (sine  space)  of  the 

° IU,IV  beam 


8.3.2  Output  Data 


The  following  data  is  output  via  the  COINCIDENCE  DETECTION  MARKS 

file: 


CDMRKS(L.M) 


This  is  the  same  file  listed  under 
Input  Data,  with  angle  data  added 
to  it. 

U I ■ 3 is  the  U coordinate 

V I ■ 4 is  the  V coordinate 
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8.4  MATHEMATICAL  RELATIONSHIPS 


Measurements  of  angular  position  are  obtained  from  monopulse  data 
for  each  marie  which  passes  coincidence  detection.  These  marks  are 
contained  In  the  Input  array  CDMARKS. 

For  each  mark  the  range  Index  k of  the  closest  range  sample  Is 
computed  from  the  range  mark  R 

( R - R ) 

K ■ Integer  < — - — — ! (8.1) 


where  Integer  {J}  means  the  Integer  closest  to  J. 

For  each  mark,  the  range-rate  index  I of  the  closest  range-rate 
sample  Is  computed  from  range-rate  mark  R 


I 


Integer 


(8.2) 


Each  mark  In  the  CDMARKS  array  contains  a pointer,  1^,  Indicating 
which  burst  (1  or  2)  produced  the  mark.  The  object's  U and  V positions 
within  the  beam  Is  then  computed  as 


4U 


dV 


Vb^  rk.i)  . 1^ 


1 

2 

1 

2 


(8.3) 


(8.4) 
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The  U and  V coordinates  of  the  beam  center,  U and  V , are  obtained 

o o 

from  a prestored  array,  based  on  the  beam  indices,  lU  and  IV. 


The  object  U and  V coordinates  are  then  computed  as 


U - U + AU 
o 


(8.5) 


V - V + AV 
o 


(8.6) 


8.5  DETAILED  DESCRIPTIONS 


8.5.1  Oblect  Beam  Position  Algorithm  Executive 
Subroutine  (ANGALG) 

8. 5. 1.1  Purpose 

The  executive  subroutine  ANGALG  controls  the  process  of  computing 
object  angular  positions. 

8. 5.1.2  Detailed  Description 

The  detailed  flow  chart  of  the  executive  subroutine  is  shown  in 
Figure  8.2.  The  computations  follow  the  sequence  described  in  Section  8.4. 
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9,0  KNOWN  OBJECT  RECOGNITION  ALGORITHM 

9.1  PURPOSE  OF  SUBSYSTEM 

The  function  of  the  Known  Object  Recognition  (KOR)  Algorithm 
is  to  identify  a target  designated  by  the  CDG  algorithm  as  one  that  is 
already  in  track  and  thus  avoid  initiation  of  a new  track  on  a previously 
designated  target.  The  track  state  may  have  resulted  from  bulk  filtering 
in  the  same  beam  or  in  a different  beam  in  an  earlier  scan.  KOR  Algorithm 
is  also  responsible  for  identifying  two  designations  of  the  same  object, 
resulting  from  the  bulk  filtering  of  two  adjacent  beams,  performed  within 
less  than  50  milliseconds  of  each  other. 

9 . 2 GENERAL  DESCRIPTION 

The  KOR  algorithm  operates  on  two  state  estimates,  one  provided  by 
the  CDG  Algorithm  and  the  other  by  the  Track  Function  (or  by  designation 
of  an  object  in  another  beam).  The  covariances  of  the  corresponding 
estimate  errors  are  also  available  to  the  KOR.  The  algorithm  consists  of 
a teat  which  decides  whether  the  two  state  estimates  (often  corresponding 
to  two  different  time  instants)  correspond  to  the  same  object  or  to  two 
different  objects.  Clearly,  this  test  should  also  use  the  covariance  of 
the  errors  in  the  designation  and  track  state  estimateion. 

The  performance  of  any  KOR  algorithm  is  specified  in  terms  of  (1)  the 
probability  of  deciding  that  the  two  state  estimates  correspond  to  the 
same  object  when  this  is  in  fact  true  and  (2)  the  probability  of  deciding 
that  the  state  estimates  correspond  to  the  same  object  when,  in  fact, 
they  correspond  to  two  different  objects.  While  the  latter  probability 
should  be  kept  very  small,  the  former  has  to  have  a reasonably  high  value 
so  that  there  is  a low  probability  of  initiating  a redundant  track. 

An  overview  flowchart  for  the  KOR  Algorithm  Is  shown  in  Figure  9.1. 
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PERFORM  X TEST  FOR  KOR  USING  \ 
A PAIR  OF  STATE  ESTIMATES,  THEIR 
ESTIMATION  TIMES  AND  ERROR 
COVARIANCES.  ONE  ELEMENT  OF  THE 
PAIR  IS  FROM  THE  COG  TRACK  FILE  AND 
THE  OTHER  IS  EITHER  FROM  THE  OBJECT 
TRACK  FILE  OR  THE  DESIGNATION  FILE. 


RETURN 


r MOST  RECENT  STATE 
ESTI.1ATES  OF  THE 
OBJECTS  IN  TRACK; 
TIME  OF  ESTIf-IATION 
AND  COVARIANCE 
^MATRICES 


(^DESIGNATION  IN  THE 
BEAMS  THAT  HAVE 
ALREADY  BEEN  | 

PROCESSED  IN  THE 
CURRENT  SCAN  OF  THE 
GROUP  OF  BEAMS 


Figure  9.1  Overview  Flowchart  of  Che  KOR  Subsystem 


An  overview  flowchart  for  the  KOR  Algorithm  is  shown  in  Figure  9.1. 


When  the  designation  sequence  in  a beam  is  completed,  the  KOR  Algorithm 
is  called.  The  input  to  KOR  consists  of  the  designations  made  up  to  that 
point,  in  the  current  scan  of  the  beam  group  and  the  state  estimates 
provided  by  the  track  function. 

Each  one  of  the  new  state  estimates  in  the  beam  being  processed  is 
paired  up  with  every  one  of  the  designations  from  the  adjacent  beams 
and  the  track  state  estimates. 

In  Figure  9.1  the  OBJECT  TRACK  file  contains  the  track  state  estimates 
along  with  estimation  covariances  and  estimation  times  of  all  of  the 
objects  that  are  being  tracked. 

The  DESIGNATION  file  contains  the  designations  (states  and  covariances) 
accepted  in  the  current  scan.  Each  designation  entry  contains  a reference 
to  the  beam  index  in  which  the  estimate  was  generated.  The  KOR  algorithm 
uses  all  of  the  state  estimates  provided  by  the  track  functior  and  only 
those  estimates  from  the  designations  of  the  current  scan  that  correspond 
to  the  beams  adjacent  to  the  one  being  processed  by  KOR.  This  is  why 
track  estimates  are  stored  in  an  .object-oriented  file  while  the  current 
scan  designations  are  stored  in  a beam-oriented  file. 

From  the  new  designations  in  the  beam  being  processed,  only  those  that 
2 

are  accepted  by  the  x test  as  new  objects  are  retained  and  added  to  the 
list  of  the  designations  accepted  in  the  current  scan.  The  rest  are 
dropped. 
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9.3  DATA  SPECIFICATION 


9.3.1  Input  Data 


The  following  data  Is  Input  from  the  CDG  Algorithm  via  the  CDG  Track 


Program 


Engineering 

Symbol 


Description 


CDGTRK  (K,M) 


DTjj(n.ni) 


Xa(1> 

x:<« 


p^(i.i) 

P^(l,2) 


P^(2,2) 

Xa(3) 

Xa(4) 


Number  of  objects  designated  In 
the  current  beam  by  the  CDG 
Algorithm  after  N transmissions 

Range-ordered  designation  list, 

Indexed  by  k. 

m * 1 Range  estimate  (m) 
m * 2 Range-rate  estimate  (m/sec) 
m > 3 Number  of  misses 

2 

m * 4 Variance  of  range  estimate  (m  ) 

m ■ 5 Covariance  of  range  and 

range-rate  estimates  (m^/sec) 

m ■ 6 Variance  of  range-rate 
estimate  (m2/sec2) 

m ■ 7 Angular  coordinate  U 

m ■ 8 Angular  coordinate  V 


The  following  data  Is  Input  via  the  DESIGNATION  file: 


desic(n;m) 


X a) 

Xg(2) 

pga.i) 

Pb(1,2) 

Pb(2.2) 


Number  of  objects  In  the  array  DESIG 

List  of  objects.  Indexed  by  N,  designated 
In  the  beams  already  processed  in  tne 
current  scan  of  the  group  of  beams. 

m > 1 Range  estimate  (m) 

m ■ 2 Range-rate  estimate  (m/sec) 

2 

0-3  Variance  of  range  estimate  (m  ) 

m ■ 4 Covariance  of  range  and 

range-rate  estimates  (m2/sec) 

m ■ 5 Variance  of  range-rate 
estimate  (m^/sec^) 


9.3.1  (Conclnued) 


Th*  following 
OBJECT  TRACK  fil«: 

NTRK 

TRACKS  (N,M) 


The  following 
GAMMA 


The  following 

rcuR 

lU(lCUR) 

IV(ICUR) 

TCDG(ICUR) 


XjO) 

m 

- 6 

Angular  coordinate  U 

XgCM 

m 

- 7 

Angular  coordinate  V 

lu 

m 

- 8 

U index  for  beam  of  designation 

IV 

m 

- 9 

V index  for  beam  of  designation 

input  fro- 

Che 

TRACK 

FUNCTION  via  the 

Number  of  objects  In  the  array  TRACKS 

List  of  Che  most  recent  data  on  Che 
objects  in  track.  Indexed  by  K 


Xg^l) 

m 

- 

1 

Range  estimate  (m) 

Xb(2) 

m 

- 

2 

Range-race  estimate  (m/sec) 

PjCl,!) 

Pb(1.2) 

m 

• 

3 

Variance  of  range  estimate  (m  ) 

m 

m 

4 

Covariance  of  range  and^ 
range-rate  estimates  (m  /sec) 

Pb(2,2) 

m 

■ 

5 

Variance  of  the  range-rate 
estimate  (m2/sec2) 

X-(3) 

m 

m 

6 

Angular  coordinate  U 

m 

m 

7 

Angular  coordinate  V 

m 

m 

8 

Time  of  estimation  (msec) 

data  is  input  via  the  ROR  PARAMETERS  file: 


Y A positive  real  number  specifying  Che 

threshold  of  the  test  performed 
for  KOR. 


data  is  input  via  the  BEAM  STATUS  file: 


lu 

Iv 

T 

COG 


Index  to  array  entries  for  current  beam 

Index  for  beam  U coordinate 

Index  for  beam  v coordinate 

Designation  time  for  current  beam 
(msec) 
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The  following  data  is  input  via  the  MEASUREMENT  PARAMETERS  file: 


SIGUV 


a ,a  Angle  measurement  standard  deviation 

U V 

for  u and  v 


9.3.2  Output  Data 

The  output  data  consists  of  an  updated  DESIGNATIONS  file  for  the 
beam  group.  Nonredundant  designations  from  the  current  beam  are  added 
to  the  file  by  the  ROR  Algorithm.  The  data  is  passed  to  the  Track 
Function  and  represents  the  principal  output  of  the  ALF. 


9.3.3  Parameter  Settings 

The  only  parameter  that  should  be  specified  is  the  threshold  y . 
Normally  the  value  of  y is  chosen  such  that  when  the  two  state  estimates 
used  in  the  KOR  test  correspond  to  two  different  objects  and  the  spacing 
of  the  two  objects  is  equal  to  or  larger  than  the  average  resolution  of 
the  designation  and  track  functions  then  the  probability  of  deciding  that 
the  two  objects  are  one  object  is  very  small.  If  this  leakage  probability 
is  assigned  a value  of  .005,  the  appropriate  y value  is  4.7. 


9.4  MATHEMATICAL  RELATIONSHIPS 


2 

The  test  for  KOR  is  a x test  of  hypothesis.  Let  x.  and  Xn 

A O 

represent  the  pair  of  state  estimate  vectors  for  which  the  test  is  to  be 

performed.  Let  P.,  t , and  P , t be  their  error  covariance  matrices  and 
A A o o 

times  of  estimation,  respectively. 

The  test  of  the  hypothesis  that  the  two  state  estimates  correspond  to 
the  same  object  consists  of  comparing  the  quadratic  form 

z A (Xb  - 

with  the  threshold  y . In  the  expression  (9.1),  we  have 

<>o  ■ ‘'b  * ‘‘H-V  ^ 

where  (t_,t  ) is  the  state  transition  matrix  which  is  uniquely  determined 
D A 

by  t_-t,  (only  constant  velocity  trajectories  are  assumed).  Therefore,  we  have 
B A 


r ^ I ^ 


L 0 . I 


(9.3) 


where  I is  the  identity  matrix.  If 
Z > Y 


(9.4) 


it  is  decided  that  the  two  state  estimates  correspond  to  two  different 
objects  and  hence  the  estimate  from  the  new  designation  is  retained. 
Otheivise,  the  two  estimates  are  assumed  to  correspond  to  the  same  object 
and  the  new  designation  la  dropped.  The  value  of  y is  a fixed  input 
parameter. 
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9.5  DETAILED  DESCRIPTIONS 


9.5.1  Executive  Subroutine 

9. 5. 1.1  Purpose 

The  EXECUTIVE  SUBROUTINE  Is  basically  responsible  for  reading  the 
appropriate  data  from  the  Input  files  and  extracting  the  pair  of  designa- 
tions which  are  passed  to  the  subroutine  CHITST.  It  also  prepares  che 
Information  that  Is  used  to  update  che  DESIGNATION  file. 

9. 5. 1.2  Detailed  Flowchart 

The  flowchart  of  KORALG  shown  In  Figure  9.2  Is  self  explanatory.  Note 
chat  since  square  beam  packing  Is  assumed  (In  u-v  space)  and  che  beams  are 
scanned  each  row  at  a time  In  che  direction  of  Increasing  v,  che  only  beams 
adjacent  Co  (IU,IV)  for  which  che  designation  sequence  Is  completed  are 
(IU-1,  IV-1),  (IU-1,IV),  (IU-1,  IV+1)  and  (lU,  IV-1) . Each  time  che  y}  test 
for  ROR  decides  Chat  che  two  designations  (one  of  them  Is  from  che  CDG 
TRAGC  file)  correspond  to  the  same  object,  that  designation  Is  deleted 
from  the  CDG  TRACK  file.  A check  la  made  to  see  whether  che  list 
has  been  diminished  to  an  empty  list.  In  which  case  control  Is  returned 
Co  che  ALF  Executive. 

When  the  CDG  TRACK  file  has  been  tested  against  all  designations  In 
adjacent  beams  and  all  OBJECT  TRACK  file  state  estimates,  che  remaining 
entries  In  che  CDG  crack  file  are  transferred  to  che  DESIGNATION  file. 

This  process  updates  che  DESIGNATION  file  for  use  on  the  next  execution 
of  che  KOR  Algorithm. 
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READ  THE  LIST  OF  DESIGNATIONS  IN 
THE  BEAM  BEING  PROCESSED. 


READ  TIME  DF  DESIGNATION  AND  THE 
BEAM  INDEX  (lU,  IV) 


IS  DESIGNATION 
COMPLETE  IN  BEAM 
(k.  M 


READ  ALL  DESIGNATIONS  FOR 
BEAM  (k,  4)  INTO  ARRAY 
BUFLST 


FOR  EACH  BEAM  (k.i)  = (IU-1  JlV-1 ) ( IU-1  .IV) ,( IU-MV+1 ) ,( lU. IV-1 ) 


OBTAIN  DESIGNATION  STATUS 
AND  DESIGNATIDN  TIME  FROM 
BEAM  STATUS  FILE  FOR  BEAM 
(k,  l) 


Figure  9.2  Known  Object  Recognition  Algorithai 
Detailed  Flowchart 
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Ftgur«  9.2  Known  Obj«cC  R«co<cnicion  Algorithm  Detailed  Flowchart 
(Continued) 
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9.5.2  Subroutine  CHITST 


9. 5. 2.1  Purpose 

2 

Subroutine  CHITST  performs  a x test  to  determine  whether  two 
designations  are  derived  from  the  same  object. 

9.5. 2. 2 Input /Outuput  Data 
Input 

^A’  Peit  of  estimated  state  vectors 

P^,  Pg  Pair  of  error  covariance  matrices 

t , t_  Pair  of  designation  times 

Ci  O 

2 

Y Threshold  for  the  x test 

Output 

ISAME  A Boolean  variable  Indicating  whether  or  not  the 

test  decided  that  designations  correspond  to  the 
same  object 

ISAME  * 1 Same  object 

ISAME  ■ 0 Two  different  objects 

9. 5. 2. 3 Detailed  Flowchart 

The  detailed  flowchart  of  CHITST  is  in  Figure  9.3.  The  subroutine 
computes  the  quadratic  form  given  by  (9.1),  (9.2)  and  (9.3)  for  the  pair 
of  designations  supplied  to  it.  The  quadratic  form  is  compared  with  y 
to  decide  whether  the  designations  are  derived  from  the  same  object. 
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f 10.0  MEMORY  REQUIREMENTS 

: ' 10.1  PURPOSE 

I 

This  section  provides  estimates  of  memory  requirements  for  the  ALF 
executive,  component  algorithms,  and  large  arrays  of  data.  These  estimates 
are  intended  only  as  initial  guidelines  for  development  programmers  when 
allocating  memory  resources . 

10.2  PROGRAM  SIZE  ESTIMATES 

Table  10.1  lists  estimated  program  size  for  the  ALF  executive 
program  and  each  of  the  component  algorithms.  The  term  "program  size" 
means  words  of  instructions  and  internal  data,  such  as  parameters,  pointers, 
small  arrays,  etc.  The  term  "program  size"  does  not  include  large  arrays 
I of  data,  such  as  radar  returns,  etc.  Large  arrays  are  described  in 

' Section  10.3. 

I Some  of  the  program  size  estimates  were  obtained  from  existing  FORTRAN 

programs  implemented  on  the  CDC  7600.  These  programs  are  marked  by  an 
asterisk  in  Table  10.1.  The  estimates  based  on  actual  programs  may  be 
excessively  large,  since  each  program  was  developed  for  use  in  a stand- 
alone mode,  i.e.,  as  a separate  execution.  Therefore,  the  development 
of  compact  code  was  not  a primary  consideration. 

Those  program  size  estimates  not  derived  from  existing  programs 
were  obtained  based  on  relative  algorithm  complexity  as  compared  to 
existing  programs. 


I 

10-1 

I 


Table  10.1 

PROGRAM  SIZE  ESTIMATES  FOR  ALF 


Program  Size 
(Decimal  Words) 


Executive  Program 

— 

1000 

Waveform  Request  Algorithm 

500 

1 

Parameter  Selection  Algorithm 

250 

it 

Marking  Algorithm 

2500 

* 

Coincidence  Detection  Algorithm 

250  1 

it 

Coherent  Double  Gating  Algorithm 

1000 

Object  Beam  Position  Algorithm 

500 

Known  Object  Recognition  Algorithm 

500 

TOTAL  ALF  SIZE 

6500 

it 

Based  on  actual  implementations 
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10.3  array  sizes 


Table  10.2  contains  size  estimates  for  the  large  arrays  of  data 
used  by  ALF.  The  file  names  and  array  names  correspond  to  the 
terminology  used  In  the  Data  Specification  subsections  for  each  algorithm. 

The  values  listed  for  Estimated  Array  Sizes  are  based  on  the  assumed 
dimensions  shown  In  the  second  column  of  Table  10.2.  These  dimensions 
are  preliminary  estimates  and  may  require  adjustments  based  on  operating 
experience  with  the  ALF  simulation. 

Some  of  the  data  files  listed  In  Table  10.2  must  be  retained  during 
the  time  between  processing  successive  transmissions  for  an  active 
beam,  l.e.,  a beam  undergoing  designation.  Data  from  other  beams  will 
normally  be  processed  during  this  time.  The  fourth  column  of  Table  10.2 
lists  the  total  number  of  versions  of  each  file  which  must  be  retained 
between  ALF  executions.  For  example,  the  BURST  1 M^RKS  file  for  an  active 
beam  must  be  saved  for  use  by  the  Coincidence  Detection  Algorithm  when  process- 
ing burst  2.  After  completion  of  burst  2 processing,  there  is  no  further 
need  to  save  the  BURST  1 MARKS  or  the  BURST  2 MARKS  data. 

Some  type  of  file  management  scheme  may  be  required  to  prevent 
the  accumulation  of  large  numbers  of  sizable,  unnecessary  data  files 
during  a multi-beam  execution  of  the  ALF. 

The  item  NBM  in  Table  10.2  represents  the  maximum  number  of  beams 
which  will  undergo  designation  processing  during  a simulated  attack. 

This  parameter  Is  highly  dependent  on  scenario  and  beamwldth,  and  Is 
not  specified  in  this  report. 
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Table  10.2 

ARRAY  SIZE  liSTIMAlXS  FOR  ALF 


File  Name 

Array  Names  and 
Dimensions 

Estimated 
Array  Sizes 

Total  Number 
of  File 
Versions  Saved 

RURST  1 VIDEO 

VIDE01(2000,9) 

2000  range  bins 

9 ranr.c-rate  bins 

18,000 

One  per  active 
beam 

BURST  2 VIDEO 

VIDE02(2000,9) 
Same  as  VIDEOl 

18,000 

None 

BURST  1 MONOPULSE 

UM0N0K2000,  9) 
Same  as  VIDEOl 
VMON01(2000,9) 
Some  as  VIDEOl 

18,000 

18,000 

One  per  active 
beam 

BURST  2 MONOPLT.se 

UM0N02  (2000,9) 
Same  as  VIDEOl 
\'MON02(2000,9) 
Same  as  VIDE02 

18,000 

18,000 

None 

BURST  1 MAR.KS 

RRDOTK800.3) 
800  marRs 
3 data  teems 
per  mark 

2,400 

One  per  active 
beam 

BURST  2 MAPJCS 

RRDOT2(800,3) 
Same  as  RRDOTl 

2,400  , 

None 

COINCIDENCE 
DETECTION  MARKS 

CDMPJCS(800,5) 

800  marks 

3 data  teems  per 
mark 

4,000 

None 

CDG  TRACK 

CDGTRK(600,3) 

600  tntctal  object 
8 data  teems  per 
object 

4,800 

One  per  active 
beam 

BEAM  DESIGNATIONS 

DESIG(50xtIBM,9) 

SO  destgnaetons  per  beam 
NBM  acctve  beams 
9 data  teems  per  design. 

450  X NBM 

One  per  desig- 
nated beam 

waveform  parameters 

BEAMS(IU,IV,4) 
lU.IV  defines  1500 
beams 

4 data  teems  per  beam 

6,000 

One 

BEAM  STATUS 

IDES(NBM),  TBURST(NBM). 
lU(NBM),  IV(NUM), 
NBPAIR(NBH) , NBUUSTCNDM) 
NBM  active  be.ams 

6 X NBM 

One 

OBJECT  TPACK 

TRACKS  (25  X NBM,  8 ) 
25  objects  per  beam 
NBM  active  beams 
8 data  ttona  per  object 

200  x'NBM 

One  [ 

UBU\M(IU,IV) 

VBILVMUU.LV) 

IU,IV  defines  1500  beams 

1500 

1500 

One 
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